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SYNTHESIS OF C"“ CARBOXYL-LABELED PHENYLALANINE! 
By G. O. HENNEBERRY, W. F. OLIVER AND B. E. BAKER 


Abstract 


C4 carboxyl-labeled phenylalanine has been synthesized by a modified Strecker 
reaction. A yield of 27.6% based on sodium carbonate (C™) has been obtained 
and the distribution of the remainder of the activity in the by-products de- 
termined. 


Introduction 


With a view to using carboxyl-labeled phenylalanine in metabolic studies 
the literature was surveyed for a feasible method of synthesizing this com- 
pound. Gurin and Delluva (2) reported a method for preparing doubly 
labeled phenylalanine. Kendall and MacKenzie (5) prepared alanine by a 
modified Strecker reaction which involved the coupling of acetaldehyde with 
ammonium chloride and sodium cyanide followed by hydrolysis. This seemed 
to be a method of making phenylalanine which could be adapted easily to the 
preparation of this compound labeled only in the carboxyl group. 


Experimental 
I. Preliminary Experiments with Nonactive Materials 


Several attempts were made to prepare phenylalanine from phenylacetalde- 
hyde by the application of the procedure of Kendall and MacKenzie, but no. 
phenylalanine was obtained. However, a small yield (6-7%) was obtained by 
following the directions given by Hickinbottom (4) for the preparation of 
phenylaminoacetic acid. In this procedure methanol was used as solvent so 
that the reactants were in a single phase. The yield was increased to 28-30% 
based on sodium cyanide by using an excess of phenylacetaldehyde freshly 
prepared by hydrolyzing 1,l-dimethoxy-2-phenylethane and by liberating 
the amino acid from the hydrochloride with silver oxide instead of with lead 
oxide. 


II. Experiments with Radioactive Carbon 

(a) Preparation of Barium Carbonate (C'*) 

The apparatus illustrated in the figure was constructed to facilitate the 
transfer of active materials and to eliminate possible contamination by atmos- 


1 Manuscript received October 12, 1950. 
Contribution from the Chemistry Department and Tracer Laboratory, Faculty of Agri- 
culture, McGill University, Macdonald College, Que. Macdonald College Journal Series No. 274. 








230 CANADIAN JOURNAL OF CHEMISTRY. VOL. 29 


pheric carbon dioxide. It consisted of a wide mouthed, 500 ml. extraction 
flask with a side arm, c, for the admission of nitrogen; three 16 mm. Pyrex 
test tubes, a, were attached to the bottom of the flask by means of 22/40 ground 
glass joints. The small desiccator, 6, was held in place by a 29/42 joint, and 
a fritted glass filter, f, was attached to a glass sleeve with DeKhotinsky cement. 











Fic. 1. Apparatus for processing active carbonates. 


The small vial containing the stock solution of sodium carbonate (C") was 
attached to the desiccator opening by means of rubber bands and 1 ml. (0.02 
mc.) transferred to each of two 16 ml. test tubes, a, with a pipette. To each 
of these tubes was then added 1 ml. of 0.4% ammonium chloride solution 
followed by 2 ml. of saturated barium hydroxide. In addition, 3 mgm. of C” 
in the form of sodium carbonate solution was added to each sample to act as 
a carrier for the barium carbonate (C™). The volume of liquid in each tube 
at this stage was approximately 7 ml. The barium carbonate precipitate was 
digested on a steam bath for one hour, cooled slowly, and separated from the 
mother liquor (Sample (A)) by means of the sintered glass filter, f. It was 
washed with water, then with acetone, and finally suspended in acetone and 
transferred to the desiccator, 6. The acetone (B) was removed by evaporation 
under reduced pressure. A yield of 141.1 mgm. of dry, very friable barium 
carbonate was obtained. 


(b) Preparation of Sodium Cyanide (C"*) 

The active carbonate was transferred to one of the reaction tubes, a, by 
means of the acetone (B) which was then evaporated and collected in a ‘“‘Dry 
Ice’’ trap. The carbonate was ground in the tube to a fine powder with a 
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snugly fitting glass pestle; 1.450 gm. of sodium azide was then added and 
mixed thoroughly with the carbonate. The mixture was decomposed at 
375°C. and then fused at 630°C. for 20 min. by a method reported previously (3) 
by two of the authors. The cyanide was distilled as hydrogen cyanide from 
the reaction mixture (Residue (C)) and trapped in about 10 ml. of water 
containing 43.2 mgm. of sodium hydroxide. A total volume of 40 ml. was 
obtained. After concentrating (Distillate (D)) the cyanide solution to approxi- 
mately 1 ml. in vacuo, the volume was made up to 25 ml. with methanol and a 
0.01 ml. aliquot (E) withdrawn for radioactive assay. The methanol solution 
was then reduced to 15 ml. by distillation under reduced pressure (Distillate 


(F)). 


(c) Preparation of Phenylalanine (C") 

The phenylacetaldehyde used in this synthesis was ted prepared by 
refluxing 10 gm. of 1,1-dimethoxy-2-phenylethane with 7 ml. of acidic 50% 
methanol for 10 min. After removing the solvents by distillation under re- 
duced pressure the crude phenylacetaldehyde was added to a solution of 
0.5575 gm. of ammonium chloride in 1.7 ml. of water and cooled to0°C. This 
was immediately followed by the methanolic solution of labeled sodium 
cyanide (15 ml.). An additional quantity of nonactive sodium cyanide 
(0.5375 gm.) dissolved in 6 ml. of 50% methanol was added very slowly to the 
reaction mixture, which was then allowed to stand at room temperature for 
18 hr. 


One milliliter of 50% sodium hydroxide solution was added to the reactants 
and the alcohol removed by distillation im vacuo (Distillate (G)). The oily 
layer was separated from the aqueous solution of inorganic salts (H) by ex- 
traction with ether. The ethereal solution was acidified with 10 ml. of 30% 
(by volume) hydrochloric acid and distilled into a solution of sodium hydroxide 
(Distillate (I)). The gummy residue was heated with three 10 ml. portions 
of 3N hydrochloric acid at the boiling point to extract the phenylalanine from 
the gummy mass. After combining the extracts the solution was evaporated’ 
to dryness under reduced pressure (Distillate (J)). It was found necessary to 
treat the crude phenylalanine with fuming hydrochloric acid at room tempera- 
ture in order to remove the last traces of gummy impurities (K). The phenyla- 
lanine solution was reduced to 2 ml. in volume; 5 ml. of ethanol was added and 
the pH adjusted to 5.96 with alcoholic ammonia. On cooling, phenylalanine 
separated out; this was then filtered off. The mother liquor was treated with 
0.1 gm. of Norite A and the volume reduced to0.5 ml. An additional amount 
of phenylalanine was obtained by adding 1.5 ml. of ethanol and cooling to 0°C. 
The combined crops of phenylalanine were recrystallized from alcohol, and a 
yield of 304.7 mgm. was obtained. The mother liquor and washings (L) were 
collected for radioactive assay. 


(d) Radioactive Assay of Samples 
All samples were monitored and only those showing activity were assayed. 
Aliquots of the samples were oxidized following the procedure of Lindenbaum 
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et al (6). However, it was found necessary to use a stronger combustion 
mixture than was recommended. It consisted of 11 gm. of chromic acid, 
67 ml. of fuming sulphuric acid (30% sulphur trioxide), and 33 ml. of syrupy 
orthophosphoric acid (sp. gr. 1.710). Each sample was treated with 17 ml. 
of this mixture. To the barium carbonate suspension from each combustion, 
1 gm. of ammonium chloride was added and the precipitate digested at 100°C. 
for one hour, then centrifuged. 


The barium carbonate plates were prepared for counting as directed by 
Armstrong and Schubert (1). The foregoing authors suggested drying the 
precipitate in a stream of air. It was found in the course of the present work 
that drying im vacuo produced more uniform precipitates. 


Two counters were used in the assay—a commercial bell-counter with a 
2.01 mgm. per sq. cm. mica window, and a bell-counter constructed in our own 
laboratory. For the latter, India ruby mica, supplied through the courtesy 
of the Mica Company of Canada, was used. It was found that this mica could 
be split under water to any desired thickness with the dexterous use of a very 
fine tungsten needle. 


The results of the assays are given in the table. The first row designates 
the samples as described in the text, the second row gives the percentage 
distribution of total activity as between the various samples. Corrections 
have been made for background, variations in counter efficiency, and self- 
absorption. Samples D, E, F contained cyanide. Accordingly these samples 
were precipitated with silver nitrate because chromic acid does not oxidize 
cyanide efficiently. The error introduced by assuming that the self-absorption 
of silver cyanide is the same as that of barium carbonate will be very small, as 
less than 2% of the total activity is involved. 








Distribution of C™ in the products of the phenylalanine synthesis 





















































Sample A|B|C|D*| E* F*| G | H| 1} j|xK KI L [Phenylalanine 
Distribution, % | 0 0. 1} 0 0 | 1.8)0.05{ 0.5] 0.2)58.91 2|0.2|10.9) 3 9) 3.7 7| 27.6 
| | 
* Samples D, E, F (cyanide). 
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THE MERCURY PHOTOSENSITIZED REACTIONS OF BENZENE 
AT HIGH TEMPERATURES! 


By E. J. Y. Scott? AND E. W. R. STEACIE 


Abstract 


An investigation has been made of the mercury photosensitized decomposition 
of benzene at high temperatures. Practically no reaction occurs at room 
temperature. At higher temperatures the main product is diphenyl although 
even at 400°C. the maximum value of diphenyl is 0.1. Ediphenyl has been found 
to be 13 kcal. mole—!. There is evidence that an activated molecule mechanism 
occurs. The secondary reactions are complex and it is not possible to arrive at a 
detailed mechanism, but the probable main reaction steps have been pointed out. 


Introduction 


The work of Bates and Taylor (1) provided the first indication that benzene 
might undergo a photosensitized reaction with mercury *P; atoms. When a 
mixture of mercury and benzene vapors was exposed to radiation from a 
mercury vapor lamp, small but definite amounts of hydrogen and methane 
were formed in addition to a tarry residue possessing an odor of diphenyl. 
Forbes and Cline (6) later produced evidence that-!P; rather than *P; atoms 
were responsible for the reaction. This was demonstrated by interposing a 
sodium nitrate filter (opaque to 1849 A but not to 2537 A radiation) between 
the cell and the lamp. Under these conditions no reaction could be detected 
although a rapid decrease in over-all pressure was observed when the filter 
was removed. 


The results described in the present paper show that although photosensit- 
ization by 2537 A produces only small conversions at room temperature, 
appreciable reaction occurs at high temperatures. 


Experimental 


Two static systems were employed depending on (a) the type of cell used, 
and (b) the position of the lamp relative to the cell. 


System I comprised the conventional static type of arrangement in which 
radiation from a mercury-neon low pressure arc was condensed by a lens onto 
the face of a cylindrical quartz cell (10 cm. long) situated in an aluminum 
block furnace. 


Considerably higher conversions were obtainable using an annular lamp-cell 
of the type shown in Fig. 1 (System II). The internal walls were constructed 
of thin quartz with a view to transmitting as much of the resonance lines as 


1 Manuscript received October 17, 1950. 
Contribution from the Division of Chemistry, National Research Council, Ottawa, Canada. 
Issued as N.R.C. No. 2329. 
2 National Research Council of Canada Postdoctorate Fellow, 1948-1950. Present address: 
Department of Chemistry, University of Notre Dame, Notre Dame, Ind. 
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Fic. 1. Lamp-reaction vessel system. 


possible. The jacket intermediate between the lamp wall and the cell served 
either as a filter compartment or as the cell proper. In order to ensure that 
the mercury vapor pressure was constant, a few drops of mercury were allowed 
to settle in or near the cell depending on whether the reaction occurred at or 
above 25°C. 


The circulatory system III was for the most part of the conventional type. 
The gases were circulated by a Puddington pump (10) through a two liter 
ballast volume, a heated mercury saturator, desaturator (at 25°C.) and cell. 
The volume of the system was 2930 cc. but when the ballast volume was not 
included in the circuit, the volume became 870 cc. The cell, consisting of a 
cylindrical quartz tube of 130 cc. capacity, 30 cm. long and 25 mm. diameter, 
was fixed in juxtaposition to a U-shaped mercury-neon low pressure arc lamp 
(11 in. long). An apparatus of this type would be expected to give a high 
intensity of resonance radiation with the minimum proportion of 1849 A. The 
latter was reduced by constructing the lamp and cell of 2 mm. thick quartz and 
by ensuring a 1 cm. air gap in the light path. The lamp reaction system was 
enclosed in an electric furnace with the electrodes projecting from the furnace. 
Temperatures were measured by a chromel-alumel thermocouple attached to 
the back face of the cell. 


Mercury cutoffs served to isolate the reaction system from the analytical, 
storage, and evacuation parts of the apparatus. Pressures could be measured 
by a manometer or by a McLeod gauge depending on their magnitude. 


Analysis 


At the conclusion of an experiment, a trap in the system was cooled with 
liquid nitrogen and circulation continued for 30 min. In this way, methane 
and hydrogen were separated from the rest of the products. Hydrogen was 
preferentially combusted on copper oxide heated to 280°C., while the remaining 
methane was analyzed, if necessary, on a Blacet—Leighton apparatus (2). In 
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experiments where a large excess of hydrogen was used, most of the hydrogen 
was removed by condensing in a silica gel trap at — 196°C. and pumping for 
30 min. 


The C.-C; fraction was then separated by a Ward—Boomer type of still (8) 
into Cs, C3, Cs, and Cx fractions. The carbon and hydrogen numbers and 
degree of unsaturation of each of these fractions could be found using the 
Blacet—Leighton apparatus. In point of fact, the amounts of these fractions 
were so small that this latter procedure was not often adopted. 


In some cases the amount of benzene consumed during the reaction was 
measured by condensing the C, + fraction into a standard volume before and 
after the experiment. 


At high temperatures diphenyl was observed as a white sublimate a few 
inches from the entrance to the furnace. This was sublimed into a trap and 
sealed off. The diphenyl was then dissolved in a standard volume of ethanol 
and estimated quantitatively by means of a Beckman Spectrophotometer 
using 1 cm. cells and 2475 A radiation. In view of the high extinction co- 
efficient of the diphenyl band at this wave length (« = 16750) small amounts 
of benzene would not invalidate these results. A sample of the white solid 
product was obtained by sublimation and found-to melt at 68-69°C. On 
mixing with pure diphenyl, no lowering of the melting point could be detected. 
Microanalysis gave: C, 93.47; H, 6.52 (Theoretical: C, 93.5; H, 6.52). 


In addition a small amount of yellow liquid was isolated, which gave an 
absorption spectrum similar to that which might be expected from a hydro- 
genated diphenyl, say Ci2 Hiz. Assuming an eof about 10,000, an approximate 
estimate of the importance of this compound could be obtained (see Table I). 


The residue in the cell would not distil at 500°C. For the purpose of obtaining 
an approximate carbon balance, this was combusted by passing in oxygen at 
550°C. for two hours and the carbon dioxide formed was collected in a trap at 
—180°C. Carbon dioxide could not be detected when a control experiment 
under identical conditions was run with no residue in the cell. Although the 
amount of carbon dioxide was only measured under one set of conditions, it 
was necessary to clean the cell by oxidizing the polymer in this fashion at the 
completion of every run. After re-evacuation, the cell was flushed with 
benzene and evacuated overnight. 


Materials 


Merck reagent benzene, stored over sodium wire, was used throughout and 
further purification was effected by bulb to bulb distillation. Hydrogen was 
taken from a commercial cylinder and purified by diffusion through a palladium 
tube. The propane, employed in the actinometric experiments, was obtained 
in cylinders from Ohio Chemical and Manufacturing Company and was also 
subjected to bulb to bulb distillation. 
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Results 


Preliminary runs were designed to verify the experimental results already 
reported in the literature. Using static system I, after 20 hr. exposure at 
25°C. with a benzene pressure of 50 mm., only 4 cu. mm. of permanent gas 
could be detected. When the temperature was raised to 400°C., after five 
hours’ illumination 48 cu. mm. were formed. Under the same conditions 
without illumination the yield of noncondensables was reduced to zero. It 
seemed safe to assume that under these conditions, only 2537 A radiation 
reached the cell. This was checked by placing a Corning filter No. 791 opaque 
to 1849 A but transmitting 50% 2537 A between the cell and the lamp. The 
rate of production of the permanent gas was approximately halved. 


Since benzene itself absorbed appreciably between 2200 and 2650 A, the 
relative importance of direct photolysis was checked by replacing the resonance 
lamp by a Hanovia S-500 lamp. This lamp emitted a reversed 2537 line but 
its intensity in the above spectral region was approximately 20 times that of 
the resonance lamp. The yield of noncondensables was found to be reduced 
by almost 90%. This afforded satisfactory evidence that the extent of photo- 
lysis was negligible compared with the photosensitized reaction. 


By using static system IJ, the experimental conditions could be altered to 
conform more closely to those used by Forbes and Cline. After 33 hr. at 25°C., 
the pressure dropped from 51 mm. to 41 mm. and a yellow residue could be 
observed on the inner wall of the annular cell. 10 cu. mm. of permanent gas 
was formed at the same time. 


In view of the above experiments it seemed certain that 2537 radiation could 
affect a photosensitized reaction at least at high temperatures. All further 
experiments were made with the circulatory system III, already described. 


Several control runs on the thermal reaction were made. At 400°C. the 
percentage of diphenyl formed thermally was found to be less than 1% of that 
formed photochemically. Similarly noncondensable gas has been measured 
up to 500°C. The effect of the thermal reaction was again found to be less than 
1% of the photochemical. 


The relative importance of the various products formed in the photosensit- 
ized reaction is shown in Table I. This is a typical circulatory run although in 
general much smaller conversions were used. In order to obtain an estimate 
of the amount of benzene consumed, the conversion in this case was large 
(62%). The reaction produces mainly diphenyl which accounts for 55% of 
the benzene decomposed. The combustible residue in the cell accounts for 
20%, hydrocarbons C,-C; for 8%, other diphenyl derivatives for 5%, while 
12% remain unaccounted for. The C2 fraction contains 50% C2He, 30% C2H,, 
and 20% CsH2 while the C; fraction comprises 85% C3Hs and 15% C3He. 


Although the amounts of hydrocarbons of low molecular weight are small, 
the significant quantities of hydrogen present in these compounds indicate that 
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TABLE I 
PRODUCTS OF THE PHOTOSENSITIZED REACTION 
Volume of system, 2930 cc. Reaction temp. 400°C. 
Volume of reaction vessel, 130 cc. Hg pressure ~2 X 10-* mm. 
Circulation rate, 175 cc./min. Exposure, 4 hr. 
Initial amount of benzene, 299 micromoles. 
CO: 
He CH, Ce C; C, Ci2H 0 Ci2Hi2 (?) from CsH¢s 
residue 
Cc. formed at N.T.P.| 0.17} 0.38} 0.39 |0.16/)0.09) — — 5.02 |- 4.18 
Micromoles formed | 7.7 | 17.1 ioe See & § 4.0 | 51.2 5.0 224 | -190:5 
Equivalent amount of 
CsHe (micromoles) 
calculated from 
carbon blance — 2.9 5.8 |3.6 |2.7 | 102.4 10.0 37.3 | -190.5 
































a considerable amount of hydrogen must be formed during the polymerization 
process. The residue is thus highly unsaturated. 


The amount of C,; to C; products may be markedly increased by the addition 
of hydrogen to the reaction as will be seen by comparison of runs 141 and 121 
in Table II. 


The quenching effects of hydrogen and benzene are of the same order of 
magnitude so that it can be assumed that in run 141 hydrogen quenches about 
90% of the radiation. These results parallel closely those of Schiff and 
Steacie (11), who found that H atoms produced by a discharge reacted with 
benzene at room temperature yielding low molecular weight hydrocarbons. 


The products found by Forbes and Cline (6) from the mercury photosensit- 
ized reaction of hydrogen—benzene mixtures at room temperature were mainly 
cyclohexadiene and hydrogenated diphenyls. It seems likely by comparison 
of runs 140, 121, and 126 that the small amounts of hydrogen and gaseous 
hydrocarbons formed at high temperature are caused by the thermal decompo- 
sition of these unsaturated compounds. 


TABLE II 


PRODUCTS OF THE PHOTOSENSITIZED REACTION OF HYDROGEN-BENZENE MIXTURES AT 400° 
Exposure: 30 min.; other conditions as in Table I 














P , Products, cc. at N.T.P. 
Run Pressure of H, 
number |benzene (mm.)| mm. He C, C; C; C, Cc; 
141 5.5 53 — 2.47 0.81 0. 65 0.19 0. 448 
121 5.0 0 0. 068 0.043 0.039 0.024 — — 
140* 7.0 42 0.090 | 0.022 0.023 0.010 — — 
126 27.0 0 0.118 0.024 0.022 0. 034 — —_ 





























* The values in this run were obtained by carrying out the reaction at 25°C. for 150 min., 
pumping off any gas formed, and heating the cell at the conclusion of the run to 400°C. 
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The investigation of the effect of pressure on the formation of diphenyl was 
limited by the vapor pressure of benzene at the trapping temperature for 
diphenyl (0°C.). Thus, only pressures between 1 and 25 mm. could be used. 
The results of these experiments are shown in Table III and Fig. 2. 


TABLE III 
EFFECT OF PRESSURE ON DIPHENYL FORMATION AT 400°C. 
Volume of system, 870 cc. Other conditions as in Table I. 





























Run Benzene Time Ci2Hi0 | Run | Benzene Time | Ci2H 10 
number | press. (mm.) | (min.) (micromoles)| number | press. (mm.) | (min.) | (micromoles) 
219 2 | 9 13.0 || 225 | 18 60 19.8 
220 5 | 30 12.0 }} 211 20 30 8.8 
222 5 | 60 27.3 || 214 20 60 | = 19.5 
224 10 60 28.3 || 229 20 90 34.0 
210 15 30 11.8 || 226 | 22.5 60 | = 17.5 
223 15 60 27.3 || 228 | 25 60 | 18.0 
227 15 2 eee Sf | 
| ' \ 
30 F— 
avs 
“a A DIPHENYL 
© Hp, 
@ CH, 
20 dn - 
‘aa 
4 x ¢; 
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Fic. 2. Effect of benzene pressure on reaction products. 


Diphenyl values taken from 60 min. runs in Table III. 
Other values from 30 min. runs in Table IV. 


Table IV and Fig. 2 show the effect of benzene pressure on other products. 


The falling off in rates at lower pressures is presumably due to incomplete 
quenching of resonance radiation. The decrease in rates of all these products 
at high pressures strongly suggests that some sort of deactivation process is 
occurring. It is of interest to note that the maximum in the curve of diphenyl 
occurs at the same pressure (10 mm.) as the corresponding curve in the ethylene 
photosensitized reaction (7), ethylene possessing a quenching effect similar to 
that of benzene. 
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TABLE IV 
EFFECT OF PRESSURE ON THE PRODUCTION OF GASEOUS PRODUCTS 
Experimental conditions as in Table I 




















| (Micromoles) 

Run Pressure of Time 

number | benzene (mm.) (min.) He CH, C2 Cs 
109 | 0.07 30 0.31 0.62 0.45 0.18 
124 0.50 15 0.27 0.98 0.36 0.18 
122 | 0.50 30 0. 67 2.14 0. 62 0.62 
123 0.50 60 0. 67 3.92 0.80 0.71 
121 5.0 30 3.06 1.92 1.74 1.07 
120 10.0 30 4.00 1.43 2.05 0. 89 
8s 27.0 15 2.85 | 0.49 0.80 — 
96 27.0 30 5.12 1.25 1.20 0.98 
126 | 27.0 30 5. 26 1.07 0.89 1.52 
97 27.0 60 11.70 2.05 3.03 3.35 
100 27.0 90 13. 80 2.94 3.35 2.54 
115 50.0 30 4.65 0.85 0.45 0.40 
114 | 94.0 30 3.12 0.45 0.54 0.36 














It will be noticed that at low pressure, methane is produced in increasing 
amounts. This phenomenon has been noted before in mercury photosensitized 
reactions (5) and it has been customary to associate this effect with atomic 
cracking reactions. The rate of methane production in the present case is 
independent of exposure time. 


Tables V and VI give the effect of temperature on the rate of disappearance 
of benzene and on the rate of formation of diphenyl. Arrhenius plots of the 
data give activation energies of 11.5 and 13 kcal. respectively. The con- 
versions were lower for the diphenyl values, and the results are therefore more 
reliable. 

TABLE V 
EFFECT OF TEMPERATURE ON THE RATE OF DISAPPEARANCE OF BENZENE 


Experimental conditions—as in Table I 
Pressure of benzene—2 mm. 











Run Temp. °C. Time Benzene reacted CieHio Kbenzene 
number (min.) (micromoles) micromoles micromoles/hr. 

185 200 240 22.1 3.30 5.5 

189 300 30 —_— 3.57 

179 300 60 18.7 5. 65 | 

177 300 60 15.6 6.89 - 17 

173 300 120 34.3 9.30 | 

172 300 240 68.7 16. 25 

181 325 30 12.5 4.80 25 

184 350 30 19.0 5.96 38 

182 375 30 22.1 6.52 44 

175 400 30 34.3 10.10 

176 400 30 50.0 9.54 | 

183 400 30 37.9 7.14 79 

187 400 30 44.2 10.09 ; 

174 400 60 87.4 14.95 | 

190 400 60 71.4 15.60 
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TABLE VI 
EFFECT OF TEMPERATURE ON THE FORMATION OF DIPHENYL 
Pressure of benzene—10 mm. 
Experimental conditions—as in Table I. 
| | ; | | | 
Diphenyl | Diphenyl | 
Run | T, °C.| Time | (micromoles | Kdipn-|| Run | T, °C.| Time | (micromoles Kdiph. 
number | (min.)} per hr.) ! number | |(min.)} perhr.) | 
171 | 200 | 60 1.5 1.5 || 163 | 300 | 360 23.7 | 4.0 
170 200 | 120 2.6 1.3 || 162 300 480 31.1 3.9 
168 200 240 3.4 0.9 | 192 350 30 6.3 12.6 
169 200 420 4.6 0.7 || 165 400 10 4.1 24.6 
191 300 60 6.7 6.7 || 166 400 | 20 8.8 | 26.4 
155 300 60 5.0 5.0 || 164 400 30 14.3 | 28.6 
160 300 120 9.4 4.7 | 188 400 30 13.0 | 26.0 
157 | 300 | 240 18.8 4.7 || 167 | 400 | 40 | 17.2 25.8 
i | 























The effect of temperature on the formation of hydrogen and gaseous hydro- 
carbons was examined at higher pressure (27 mm.) and temperatures (350 to 
500°C.). In this way the amounts of gas obtainable were increased while 
maintaining small percentage conversions. 
from these results are unexpectedly high: Ey, = 24.5 and Ecy, = 23.5 
kcal. per mole but in view of the relatively small amounts of these products 
and the complexity of the processes forming them, they are probably not of 
much significance. 


Pressure of benzene—27 


TABLE VII 


EFFECT OF TEMPERATURE ON THE PRODUCTION OF HYDROGEN AND HYDROCARBONS 


mm. 
Other conditions—as in Table I. 


The activation energies obtained 
































| Micromoles 
Run_ | Temp., °C.| Exposure | 
number | time, min. | H2 Ch, | C2 C; Cy 
91 350 120 0.38 O11 | 0:05 —_ | — 
101 350 180 0.44 0.09 | 0.10 0.10 | — 
88 400 15 0.29 0.05 0.08 — a 
87 400 30 0.59 0.10 0. 22 — 
96 400 30 0.51 0.13 0.12 020 | 0:10 
97 400 60 re ri 0.21 | 0.30 0.35 | 0.14 
100 400 90 1.38 0.29 | 0.35 0.25 | 0.04 
90 450 15 1.03 0.21 0.12 0.08 — 
89 450 30 1.85 0. 29 0. 29 0.25 | — 
99 450 60 3.29 0. 64 0.70 0.43 | 0.08 
92 500 5 0.97 0.13 — — — 
98 500 10 2.14 0.30 0. 24 0.10 | 0. 003 
94 500 15 3.73 0.32 | 0.06 0.05 — 
103 | 500 20 3. 82 0.74 | 0.83 0.34 | 0.09 
i 





Quantum Efficiency 


Since most of the experiments were carried out in the region of 400°C. 


actinometric experiments were made at that temperature. 


The quantum 


yield of hydrogen formation in the mercury photosensitized reaction of propane 
has been found by Bywater and Steacie (4) to be 1.4 at 140 mm. pressure and 
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400°C. Using their value, the intensity of radiation was found to be 2.8 X 
10~ Einstein per hr. corresponding to a ¢gipheny: Of 0.1 at 10 mm. pressure. 
The quantum yield of benzene disappearance is therefore of the order of 0.2. 


These values are only to be considered approximate but they show that, 
even under the optimum conditions, the reaction has a low quantum yield. 
This supports the view that some kind of deactivating process occurs. 


Discussion 


The diverse nature of the products indicates that the mercury photosensitized 
reaction of benzene is extremely complex. However, many of the products, 
such as hydrogen, C,; to Cy, hydrocarbons, and probably the cell residue, may 
be presumed to be the result of secondary reactions at the high temperatures 
employed. If diphenyl is considered to be the most significant product, .as 
seems to be the case at least at low pressures, then two significant facts must 
be explained; first, the lack of reaction at room temperature and second, the 
activation energy of approximately 13 kcal. for diphenyl formation and benzene 

disappearance. 


The effect of pressure on the reaction is quite definite. Increase of pressure 
causes a decrease in yield of all the reaction products provided quenching is 
complete. Apparently in the initial step an activated molecule is produced, 


Hg(*?P1) + CeHe = Hg(*So) + CoHe™, (1) 
which is deactivated at higher pressures by (2). 
C.H.e* + M=CeHe + M. . (2) 


Since dipheny] is the main reaction product it is probably formed by an initial 

rupture of the C—-H bond followed by mutual combination of the phenyl 
radicals generated. 

C.H,* = C.H; + H (3) 

2C.Hs = CuHw. (4) 


The above steps account for the effect of pressure on the formation of 
dipheny] but are not sufficient to explain the lack of reaction at room tempera- 
ture. 


It seems most unlikely that Reactions (1), (2), or (3) should have an apprec- 
iable activation energy. It, therefore, seems certain that the lack of reaction at 
room temperature must be due to a re-formation of benzene from C,H; and H. 
It is known from a variety of evidence that H-atoms react rapidly with 
benzene (3, 6, 9, 11). The mechanism of re-formation of benzene might 
therefore be the reverse of Reaction (3), or 


H + CeHe = CoH (5) 
C.H; + CeHs = 2CeHe. (6) 


This is, in effect, equivalent to asserting that benzene is an unusually efficient 
third body for phenyl! radical and hydrogen atom recombination. 
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In order to account for the activation energy of benzene disappearance and 
diphenyl formation it seems essential to assume that alternative radical re- 
actions occur which compete with the recombination steps at high tempera- 
tures. The most likely reaction of this sort seems to be 


H + C.Hs = CoH; + H2 (7) 


which would be expected to have an activation energy of the order of magnitude 
of 13 kcal. 


Another alternative reaction which cannot be excluded is an attack of a 
phenyl radical on a benzene molecule by (8). 
C,H; + CeHe = Cy2H i + H . (8) 


This reaction violates the Rice—Teller principle, but seems to occur readily in 
solution (12). 


The above mechanism is by no means complete or established. It seems 
plausible, however, and there is no point in further speculation on the basis of 
existing data. 
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SOUND VELOCITY AND SOUND ABSORPTION IN 
THE CRITICAL TEMPERATURE REGION! 


By W. G. SCHNEIDER 


Abstract 


The velocity and absorption of ultrasound (600 kc.) has been measured 
throughout the critical temperature region of sulphur hexafluoride. Measure- 
ments were carried out for the coexisting liquid phase and vapor phase below T-;, 
and for the supercritical gas, and simultaneously, observations of the meniscus 
behavior in the neighborhood of 7, were made. The sound velocity for both 
liquid and vapor phases below 7, decreased with increasing temperature and 
became equal at 7;, the velocity at this point being 121.5 m. per sec. In the 
temperature range from 0. 6° below T, to T, the velocity in the vapor was greater 
than that in the liquid. A very high absorption of sound was observed, having 
a maximum value at 7, and extending over a temperature range of approxi- 
mately 1°. In the temperature range from T; to 0.6° below 7;, the absorption 
in the liquid phase was greater than that in the vapor. 


Introduction 


An anomalous and very high absorption of ultrasound in the critical tempera- 
ture region of sulphur hexafluoride has been previously reported (2). A more 
detailed series of measurements for both the coexisting liquid and vapor phase 
has now been completed. In these measurements the sound velocity and 
sound absorption were measured alternately in the liquid phase and vapor 
phase as the system was heated from a temperature a few degrees below the 
critical to a temperature several degrees above the critical temperature. The 
measurements were made with an acoustic interferometer built into a heavy- 
walled glass tube permitting observations of the meniscus and of the position 
of the quartz crystals of the interferometer. Sulphur hexafluoride was used 
in these experiments both because its molecules are nonpolar and spherically 
symmetric and because it has very convenient critical properties (7, = 
45.55°C., P. = 36 atm.). 


Experimental 


Description of the Acoustic Interferometer 

The acoustic interferometer, schematically shown in Fig. 1, was of the 
double crystal type, and was in many respects similar to the interferometer 
previously described (3). The acoustic chamber consisted of a heavy-walled 
cylindrical glass tube A, 3 in. long and 15/16 in. internal diameter. Pressure 
tight seals were obtained on the squared ends of the tube by means of lead 
washers compressed on the glass cylinder from both ends by plunger C. This 
plunger also houses the keyed piston P, on which the quartz crystal is mounted, 
and the metric screw drive. In order to assure perfect alignment of the 

1 Manuscript received, October 17, 1950. 
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Fic. 1. Schematic diagram of the acoustic interferometer. 


crystals, the plunger C is machined to a snug fit in the outer cylinder B. The 
transmitting crystal T is mounted on the lower piston and the receiving crystal 
on the upper piston. The crystal leads were brought out through capillary 
tubes and through the pressures seals S. To the upper pressure seal a steel 
capillary tube was joined for evacuating and filling the interferometer. 


By means of a gear arrangement both crystals could be moved into the 
liquid phase or into the vapor phase. The receiver crystal could also be moved 
at a slow and uniform rate by the motor M, simply by inserting a tapered pin 
into the gear drive. The auxiliary apparatus and the procedure for the 
velocity of sound measurements were exactly as described previously (3). By 
means of a window slotted into the chamber B, the whole length of the acoustic 
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chamber was visible and the position of the crystals as well as the position and 
behavior of the meniscus could be observed at all times. 


All measurements were made at a frequency of 600 kc. 


Procedure 


Before assembly all parts of the interferometer were thoroughly cleaned. 
After mounting in a thermostat, capable of maintaining the temperature 
constant to within 0.002°C., the interferometer was thoroughly evacuated 
and flushed with pure sulphur hexafluoride. To fill the interferometer to a 
density close to the critical density, the thermostat was held at a temperature 
approximately 0.5° below the critical temperature, and gas was admitted 
slowly from a storage bomb and allowed to condense in the interferometer 
until the liquid level was exactly halfway up the glass tube A. That the 
filling obtained in this manner was close to the critical density was evidenced 
by the fact that, on slowly raising the temperature, the meniscus disappeared 
without perceptible motion up or down in the tube. All measurements re- 
ported in this paper were made on a single filling of the interferometer. It is 
perhaps also worth emphasizing that the system was not stirred mechanically 
‘at any time (apart from the gentle mixing caused by movement of the crystals), 
and the critical temperature, 7., was taken as that temperature at which the 
meniscus was observed to disappear. For the relationships between the 
critical temperature so defined, and the critical temperature which has been 
defined for a stirred system, see reference (1). 


Readings of the sound velocity were obtained in the following manner. 
Both crystals were moved into the phase to be measured and positioned about 
1 cm. apart. The thermostat was adjusted to the desired temperature and 
the system was then allowed to thermostat for 12 hr., after which a reading 
was taken by moving the receiver crystal, by means of motor M, towards or 
away from the transmitter crystal and recording the resonance peaks. The 
12 hr. thermostating period was found to be necessary to obtain reproducible 
results; readings taken at two hour and also four hour intervals were found to 
be erratic, but after 12 hr. an apparently stable equilibrium had been reached, 
the measurements remaining unaltered 24 or 36 hr. later. For the vapor the 
velocity readings were taken in a region approximately 1 cm. below the top of 
the interferometer, and for the liquid approximately 1 cm. above the bottom 
(or approximately 3 cm. below the meniscus level). 


Results 
Sound Velocity 
The velocity of sound measurements through the critical temperature region 
are summarized in Table I and plotted as a function of temperature in Figs. 2 
and 3 for the vapor and liquid phases respectively. To obtain the vapor (or 
liquid) data, readings were commenced several degrees below the critical 
temperature and the temperature was raised in steps to several degrees above 
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the critical temperature. 


As stated earlier, each reading was taken only 


after thermostating at a given temperature for 12 hr. Then to check the 
reversibility of the system, after the highest temperature was reached, the 
temperature was lowered in steps back through the critical temperature to 
the starting temperature, as shown by the black circles on the velocity curve 
of Fig. 2. A similar reversibility had previously been observed for the liquid 


phase (2). 


TABLE I 


SUMMARY OF SOUND VELOCITY AND SOUND ABSORPTION MEASUREMENTS 














Sound | 
Temp., °C. velocity, | Absorption | Remarks* 
| m./sec. | coefficient | 
Vapor phase 
41.917 151.20 0.050 T 
42.003 150.7 0.063 L 
42.732 148.8 0.055 L 
42.742 148.8 0.048 T 
43. 605 146.4 0. 056 T 
43.757 145.5 0.058 L 
44. 280 143.5 0.074 T 
44. 626 —_ 0.097 1 
44. 681 141.4 0.099 L 
44.952 — 0.111 + 
45.189 136.8 0. 132 L 
45. 388 133.7 0. 169 L 
45. 394 133.2 0. 161 T 
45.478 131.2 0.195 T 
45. 489 130.7 0.198 k 
45. 502 130.7 0. 208 
45.553 128.2 0.271 t Faint opalescence 
45.572 128.2 0. 292 | Opalescence more intense 
45.617 125.2 0.374 Tt Opalescence more intense 
45. 647 124.9 0. 410 | Heavy opal.—meniscus visib. 
45. 669 121.6 0.518 T f Very heavy brownish opal. 
Meniscus no longer visible 

45. 688 122.2 0. 528 { Meniscus no longer visible 
45.724 122.1 0. 454 T Opalescence less intense 
45. 766 123.4 0.410 | Opalescence less intense 
45.773 123.3 0. 382 Tt Opalescence fainter 
45. 839 123.9 0. 330 t Faint opalescence 
45. 861 123.7 0.328 { Faint opalescence 
45. 936 124.7 0. 266 T Very faint opalescence ' 
46. 084 127.0 0.199 Tt Very faint opalescence 
46.178 127.9 0. 184 | Very faint opalescence 
46. 329 129.9 0. 143 T 
46. 340 —_ 0. 157 E 
46. 854 134. 1 0.065 + 
47.438 137.4 0.043 J 
47.924 139.9 0.044 tT 
48.479 142.0 0. 041 
49.013 144.4 0.041 x 











* Indicates measurements were taken during series with rising temperature. 
| Indicates measurements were taken during series with decreasing temperature. 
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TABLE I (Cont'd.) 











| Sound 
Temp., °C. | velocity, Absorption Remarks* 
| m./sec. coefficient 
Liquid phase 
42.003 194.8 0. 060 t 
43.035 179.0 0. 064 7 
43. 882 163.6 0.072 Tt 
44. 698 146.7 0. 089 Tt 
45. 152 134.7 0.173 s 
45.379 129.2 0. 260 tT 
45. 487 125.5 0.360 t No visible opalescence 
45. 535 124.3 0.395 t Very faint opalescence 
45.590 122.5 0. 456 t Faint opalescence 
45.639 121.4 0. 546 t Heavy opalescence—meniscus distinguishable 
45.691 121.8 0.540 tT Very heavy opalescence—no meniscus 
45.739 122.6 0. 466 Tt Heavy opalescence 
45.785 123.3 0. 404 t+ Opalescence less intense 
45. 847 124.2 0.342 t Opalescence fainter 
45.930 125.2 0.280 T Opalescence faint 
46.053 126.4 0.221 Tt Opalescence very faint 
46. 239 128.4 0.177 Tt Opalescence very faint 
46. 687 133.1 0.120 7 
47. 482 137.8 0.049 ‘J 
48. 531 142.5 0.048 t 
49.149 145.3 0.049 tT 














* + Indicates measurements were taken during series with rising temperature. 
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Fic. 3. Plot of the sound velocity and sound absorption in the liquid phase vs. temperature. 


For comparison the velocity curves of Figs. 2 and 3, which were measured 
separately for the vapor and liquid phases, are shown plotted together in Fig. 4. 
In this plot the full curve represents the ‘‘liquid’’ curve, the individual points 
on the curve having been omitted for the sake of clarity. The circles and the 
dotted curve represent the vapor measurements. It should be pointed out 
that the branch of the velocity curve above T, for the liquid curve shown in 
Figs. 3 and 4 is actually for a supercritical gas, but differs from the corresponding 
branch (i.e., above 7.) of Fig. 2 in that for the “‘liquid’’ curve the measure- 
ments were made throughout at a point approximately 1 cm. above the 
bottom of the interferometer; for the vapor curve the corresponding readings 
were taken at the top of the interferometer. 
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Fic. 4. Superposition of sound velocity data for the liquid qnd the vapor phase; the full curve 
is for the liquid phase, the circles and dotted curve for the vapor phase. 
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As shown in Fig. 4, the velocity curves for liquid and vapor cross at a 
temperature approximately 0.6° below T,, and from this temperature up to 
the critical temperature the sound velocity is higher in the vapor than in the 
liquid. This difference, which is at present inexplicable, is well outside the 
limits of experimental error and must be assumed to be real; the excellent 
superposition of the two curves above T, lends further support to this con- 
clusion. 


Both the liquid and vapor curves exhibit a pronounced minimum in the 
velocity at T., the temperature of meniscus disappearance, which was observed 
to be 45.65° + 0.01°C. The minimum in the vapor curve occurs at 45.669°C., 
and that for the liquid curve at 45.639°C. However, since in this region the 
temperature was varied in steps of approximately 0.05°, the minima of the 
two curves can, within this limit, be considered coincident. The sound 
velocity at this point is 121.5 m. per sec. 


The sound velocity can be expressed by the relation: 


oP , 
= =o (fo aw ’ (1) 
Gi, 


where c = sound velocity, 
v = volume, 
P = pressure, 
andy = Cp/¢». 


; ‘ : oP al 
At first sight it would appear that since (*) tends to zero at the critical 
wv /T : 
point, the sound velocity should tend to zero. However y simultaneously 
tends to infinity and hence ¢ can remain finite. The combination of the two 
effects gives rise to the characteristic minimum at the critical point. 


Sound Absorption 

The sound absorption was measured at each temperature simultaneously 
with the velocity measurements, and was observed to become very large in 
the immediate neighborhood of the critical point. In Table I, and Figs. 2 
and 3, the absorption coefficient is shown as a function of temperature. For 
the present purpose we have arbitrarily defined the absorption coefficient a by 
the relation: 

A = Age ™, (2) 
where Ap is the initial amplitude of the sound wave, A the amplitude at 
distance x, and X is the wave length. The absorption coefficient a is hence 
the absorption coefficient per wave length, and was measured directly from 
the recorded resonance peaks. An example of a recording of the resonance 
peaks obtained at a temperature close to T, is shown in Fig. 5. To measure a, 
a smooth curve was drawn through the crests of the resonance peaks, and a 
base line was drawn through the troughs; the amplitudes A and Ao were then 
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Fic. 5. Reproduction of a typical recorder tracing obtained at a temperature close to T, showing 
a rapid falling-off in amplitude of the resonance peak. 


read off relative to the base line. Thus if Ao is the amplitude of the first peak 
obtained in this manner, and Ajo the amplitude of the 20th peak (since the 
peak separation is \/2), we have from Equation (2): 

ro” 

a = —Iln —.: 

10 = Ao 
The a@ obtained in this way is not the true absorption coefficient, and, since a 
double crystal interferometer was used in this work, it is not clear what relation 
it bears to the absorption coefficient as ordinarily defined. In the present 
work, however, we are concerned only with relative values of a as a function 
of temperature, and an exact analysis is not required. 


(3) 


As shown by the curves in Figs. 2 and 3, the absorption coefficient increases 
rapidly just below the critical temperature and attains a maximum value at 
T,, the point at which the meniscus disappears and the sound velocity attains 
a minimum value. Above this temperature the absorption again diminishes 
to alow value. The absorption measurements could be accurately reproduced 
with either a stepwise increase or decrease in the temperature (and with a 
12 hr. period of equilibration before each measurement) (see Fig. 2), and hence 
the resulting absorption curve is to be regarded as a characteristic property of 
a system in the critical temperature region. It is noteworthy, in this con- 
nection, that the temperature range over which high sound absorption occurs 
is approximately 1°, and extends over a range both below and above the 
critical point. 


Fig. 6 shows a superposition of the absorption curves which were measured 
separately for the liquid and vapor phases. As was the case with the sound 
velocity curve, the branches of the absorption curve above T, superimpose to 





SCHNEIDER: SOUND VELOCITY AND SOUND ABSORPTION 251 


9 
o 





LIQUID PHASE 


| VAPOUR aie, 


a 


ore ero 


2 
rs 
| 






SOUND ABSORPTION COEFFICIENT 
ny 
| 








9° 
° 
T 


l ; l ; l P l 
42 44 46 48 50 


TEMPERATURE, °C. 


_ Fic. 6. Superposition of the sound absorption data for the liquid and the vapor phase; the full 
curve ts for the liquid phase, the circles and the dotted curve for the vapor phase. 








within the experimental error. Over a temperature interval extending to 
approximately 0.6° below T., the absorption curve for the liquid is considera- 
bly higher. We noted earlier that in this temperature region the sound 
velocity curves for the two phases crossed, the vapor exhibiting the higher 
velocity. 


To explain the high sound absorption in the critical region two alternative 
explanations suggest themselves. The first of these is based on a simple 
diffraction of the sound beam, the second on the behavior of a relaxation 
phenomenon which is characteristic of the system in this region. 


To consider the first of these two alternatives we may regard the sound bear 
as undergoing a ‘“‘scattering’’, similar to the scattering of light in the critical 
region. The latter is due to the characteristic critical opalescence, and it is 
noteworthy that the opalescence is most intense in the temperature region 
where the absorption curve exhibits a maximum. (See Table I.) In order 
to obtain scattering, the scattering “‘centers’’ must be of the order of magnitude 
of the wave length, which in our experiments near T, was approximately 
0.2 mm. To realize scattering centers of this magnitude, one must assume 
the existence of fairly large droplets in the vapor phase, and gas bubbles in 
the liquid phase. Under the influence of a gravitational field it is difficult to 
see why such a system should be stable over long periods of time. 


According to our second model we postulate the existence of a configurational 
relaxation effect. If the spatial configuration of the system is different before 
and after passage of the sound wave, an entropy change takes place which 
leads to an excess absorption of sound energy. Alternatively, if molecular 
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clustering is assumed, the clusters may be partially disrupted during the rare- 
faction phase of the sound wave, and during the ensuing compression do not 
have sufficient time to re-form, or else they are re-formed with a different 
configuration. In other words the “‘structural’’ changes taking place in the 
system follow a time dependence which causes them to be out of phase with 
the sound field. 
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COMPOSITIONAL HETEROGENEITY OF 
BUTADIENE-ACRYLONITRILE COPOLYMERS 
PREPARED IN EMULSION AT 5°C.! 


By W. H. EMBREE, J. M. MITCHELL, AND H. LEVERNE WILLIAMS 


Abstract 


The copolymerization of butadiene and acrylonitrile is very similar to the 
copolymerization of butadiene and styrene. Polymers predominantly butadiene 
may be studied by conventional solution techniques but the study of polymers 
rich in acrylonitrile requires improved solvents for these materials. Polymer- 
ization rates are greatest for monomer ratios approximating equal proportions. 
The mercaptan modifier disappears much more slowly than in the butadiene— 
styrene system, the regulating index approximating unity. The number aver- 
age molecular weights calculated from the mercaptan disappearance curves 
indicate uniform polymer molecular weights to relatively high conversions after 
which there is a decrease. The viscosity data indicate a rise in viscosity with 
conversion, which effect is overcome for charges rich in acrylonitrile by the 
lessening of branching, the more rapid disappearance of mercaptan at high con- 
version, and the tendency of polymers containing over 50% acrylonitrile to show 
very low dilute solution viscosities in the solvents tested. Vi iscosity molecular 
weights have been calculated and estimates of the molecular weight distribution 
made. These distributions appear to be quite narrow and the usual broadening 
at higher conversions is prevented by the increased modifier consumption and 
increased vinyl content of the polymer prepared with 50 parts acrylonitrile in 
the charge. The bound acrylonitrile has been determined at various con- 
versions and the reactivity ratios have been found to be 7; = 0.28 and rz = 0.02 
for emulsions and 7; = 0. 18 and rz = 0.03 for oil phase portion only. Q is 0.74 
and e is 1.47 as calculated by the Alfrey—Price equations. 


Introduction 


Many acrylonitrile polymers and copolymers have been studied but the 
data published on the emulsion copolymerization of butadiene and acrylonitrile 
are scarce. The most directly useful results are those of Wall (23, 24) and 
of Semon (19, 20). Some data obtained during the copolymerization of other 
monomers with acrylonitrile (4, 5, 6, 9, 15) may be used to calculate the 
approximate behavior of the butadiene—acrylonitrile system. The purpose of 


this work is to co-ordinate the data and to fill in the gaps in the published 
information. 


The butadiene-acrylonitrile system is complicated by two factors. The 
polymer is insoluble in butadiene which may make the oil phase portion of 
the system heterogeneous. Acrylonitrile is appreciably soluble in water and 
is extremely reactive towards many common materials such as mercaptan, 
peroxides, amines, etc. Thus the distribution of acrylonitrile between the oil 


and water phases and its interaction with polymerization recipe ingredients 
must be considered. 


1 
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Two series of experiments were conducted some years apart. In the first a 
recipe similar to one used for butadiene-styrene studies was used (see (13)). 
When interest became centered around low temperature polymerizations to 
yield better and more uniform butadiene-styrene copolymers, the activated 
recipes (see (12, 22)) were applied successfully to the butadiene—acrylonitrile 
system. 


Techniques 


The polymerization technique was that generally used in the rubber indus- 
try (3, 12, 13, 16, 22). Modifier disappearance curves and intrinsic vistex 
viscosities were determined conventionally (16). Bound acrylonitrile was 
estimated by Kjeldahl nitrogen determinations. The fractionation curves 
were prepared from the usual solvent—nonsolvent precipitation procedure (see 
(14)). Gel (3, 13) was considered as all insoluble material in the particular 
solvent, and the swelling index (3, 13) as the ratio of the weight of swollen gel 
to weight of polymer in the gel. The dilute solution viscosity (3, 13) was 
measured on the soluble portion of the dry polymer. Monomer charge ratios 
are weight to weight, which with butadiene and acrylonitrile are very nearly 
mole to mole. 


Experimental 


Dilute Solution Viscosity and Gel of Polymers Prepared at 30°C. 
The effect of degree of conversion on the solution properties of polymers is 
shown in Table I. The gel and swelling index remain relatively constant but 


TABLE I 


EFFECT OF CONVERSION ON THE GEL, SWELLING INDEX, AND DILUTE SOLUTION VISCOSITY 








Swelling index Dilute solution viscosity 





Time (hr.) % Conversion % Gel 
in chlorobenzene solution 





71.5/28.5 charge ratio butadiene to acrylonitrile 





4 37.1 49 65 2.88 
6 46.2 46 50 2.74 
8 57.3 44 82 2.64 
11 68. 2 70 67 3.49 
15 80.1 46 53 1.92 
20 86.7 41 52 0.84 





60/40 charge ratio butadiene to acrylonitrile 





3 22 12 19 1.95 
5 56.6 25 41 2. 84 
7 i eos 49 49 1.73 
10 93.7 72 38 0.91 
14 96.2 94 31 0. 86 
19 96.8 94 26 0. 82 
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the dilute solution viscosity of the soluble portion rises to a maximum and then 
decreases. Somewhat similar results are obtained for the dilute solution 
viscosity of a 60-40 butadiene—acrylonitrile charge. There tends to be an 
increase in the amount of gel with conversion more in keeping with the lower 
dilute solution viscosity. The latter may also be due to the higher acrylonitrile 
content of the polymer or reduced effectiveness of the modifier which was 
unchanged from the previous experiment, i.e., 0.45 part of dodecyl mercaptan 
per 100 parts monomers. 


Fractional Precipitation of Polymers Prepared at 30°C. 

Three of the polymers were submitted to fractional precipitation from 
chlorobenzene—methanol solvent—nonsolvent system. The data for the samples 
are shown in Table II. The gel content of the second sample appears high 











TABLE II 
SAMPLES FOR FRACTIONATION 
Sample Charge ratio % Conversion % Gel Dilute solution viscosity 
. T-551-1* 75/25 74.5 — i 
AM-4, 16-5-4 71.5/28.5 68.2 70 3.49 
AM-4, 21-6-3 60/40 ri ie J 49 1.73 














* A sample prepared in the Pilot Plant 5-gallon glass lined autoclave. 


in view of the high dilute solution viscosity. The integral and differential 
distribution curves are in Fig. 1. T-—551-—1 starts to precipitate with 45 volume 
% methanol content. There is a peak at 45.6 volume % and most of the 
polymer has precipitated when 46 volume % has been reached. AM-4, 16-5-4 
begins to precipitate with 44 volume % methanol, reaches a slightly higher 
peak at 45 volume % and is largely precipitated with 46 volume %. AM-4,21- 
6-3 being of lower viscosity does not start to precipitate until 46 volume % 
methanol has been reached, a short peak is formed at 47 volume % and the 
polymer has largely precipitated by 48.5 volume % methanol. Thus the 
behavior of the polymers was qualitatively similar to that observed in the 
butadiene-styrene system but the precipitation becomes evident only at high 
methanol contents of the mixed solvent-nonsolvent, and the molecular weight 
distribution appears very narrow. 


A thorough study of sol-gel equilibrium, swelling index, viscosity, and 
fractionation curves is required since considerable evidence in the above work 
was found for strong dependence of some of these variables on solvent to 
polymer ratio. 


Effect of Charge Ratio on the Rate of Conversion at 5°C. 

The effect of variation of the butadiene—acrylonitrile charge ratio on the 
rate of conversion was studied. The data are in Table III and show clearly 
that 63/37 and 50/50 charge ratios polymerize at the fastest rates. The cause 
for this optimal rate of conversion will be discussed later. 
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Fic. 1. Fraction of polymer in solution (F) as a function of volume per cent of methanol (@) in 
chlorobenzene-methanol mixture and differential curves (--dF/d) for sample A M-4, 16-5-4 (e), 
T-551-1 (@) and A M-4, 21-6-3 (a). 





TABLE III 
EFFECT OF CHARGE RATIO ON THE RATE OF CONVERSION AT 5°C. 














Charge ratio—butadiene/acrylonitrile | % conversion in 17 hr. 
75/25 75.8 
63/37 96.7 
50/50 93.8 
25/75 | 79.0 
0/100 0.0 





Modifier Disappearance Curves at 5°C. 


Some measure of the molecular weight heterogeneity of the polymer may be 
obtained by measuring the modifier disappearance curves. Titration of the 
mercaptan was impossible in the normal charge; even without polymerization 
there appeared to be no residual mercaptan. This was caused by a reaction 
between the acrylonitrile and the mercaptan (1), whereby the mercaptan was 
no longer titratable although it still functioned as a modifier or chain transfer 
agent. This reaction is: 


RSH + CH: = CHCN — RSCH:2 — CH2:CN 


EMBREE ET AL.: BUTADIENE-ACRYLONITRILE COPOLYMERS 257 


However, this difficulty could be overcome by adding the mercaptan, 0.7 
part MTM* per 100 parts monomer, after the acrylonitrile had become diluted 
with the butadiene and distributed between the two phases. The most con- 
venient method was to add it by syringe and needle to the completed charge 
optionally in the form of an alcoholic solution. This technique proved ade- 
quate to obtain the modifier disappearance data shown in Fig. 2. 
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Fic. 2. Modifier disappearance data for mixed tertiary mercaptan blend MTM for 90/10 
butadiene/acrylonitrile charge ratio (ms), 66/34 charge ratio (e), and 50/50 charge ratio (a). 





It is evident that the rate of disappearance is considerably slower than in 
butadiene-styrene system (3, 13) and that the rate diminishes with increasing 
acrylonitrile in the charge. Thus must be interpreted tentatively in terms 
of a low reactivity of the acrylonitrylyl radical for mercaptan hydrogen. A 
more comprehensive study of the action of mercaptans as modifiers in this 
system is to be published (7). 


The regulating index, defined as 
—dinR 
es 
dP 


where R is the per cent residual mercaptan at fractional conversion P, was 
calculated for each charge ratio. The results were 1.38 for 90/10, 0.83 for 
66/34, and 0.71 for 50/50 charge ratio. These offer a means of numerical 
comparison with similar calculations in the GR-S and Buna S-3 system (3, 13). 
There was no waste factor evident in any of the results. 


* Phillips Petroleum Co. blend of 60% tertiary C-12, 20% tertiary C-14, and 20% tertiary 
C-16 mercaptans. 
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Number Average Molecular Weights 

Assuming an average molecular weight for the mercaptan of 220 and that 
one mercaptan molecule is consumed per molecule of polymer formed, the 
number average molecular weights based on the data in Fig. 2 were calculated. 
The results are in Table IV. At low conversions where branching is small 
and the termination by other than mercaptan is minimal, the results should 
be reasonably representative. 


TABLE IV 


NUMBER AVERAGE MOLECULAR WEIGHTS CALCULATED FROM MODIFIER DISAPPEARANCE DATA 








90/10 charge ratio 66/34 charge ratio | 50/50 charge ratio 

% Conv. | | 

= Mercaptan | Molecular | Mercaptan Molecular Mercaptan | Molecular 
consumed, %| weight calc. |consumed, %| weight calc. jconsumed, %| weight calc. 











10 14 | 2.2 * 10 | 8 | 3.9 X 10 | 8 | 3.9 X 104 
20 25 | - 2.8 | 14 | 4.5 | 14 | 4.5 
30 35 2.7 21 | 4.5 | 21 4.5 
40 47 2.7 27 4.7 27 ; Be 
50 61 2.6 33 4.8 37 | 4.8 
60 722 286| (2.6 44 | 46 | 41. | 46 
70 82 2s) CC 54 4.1 54 | 4.1 
80 | 94 | 2.7 74 3.4 | 71 3.5 
90 | 98 2.9 — | — 84 3.4 





These calculations show that the molecular weights of the 90/10 charge ratio 
polymers are relatively uniform at all conversions and that the molecular weights 
of 66/34 and 50/50 charge ratio polymers are relatively constant to high con- 
versions at which time the new polymers are of lower molecular weights so that 
the average is lower. This is shown more clearly by calculating the increment 
number average molecular weights, Table V, in the same manner. The 
increment molecular weights for the 90/10 charge ratio remain constant but 
those for the 66/34 and 50/50 charge ratios decrease above 60% conversion 
after passing through a broad maximum. 


TABLE V 


CALCULATED INCREMENT NUMBER AVERAGE MOLECULAR WEIGHTS FROM 
MODIFIER DISAPPEARANCE DATA 








Increment number average molecular weight 





Conversion range 





| 90/10 charge ratio 66/34 charge ratio | 50/50 charge ratio 
0-10% 2.2 X 10¢ 3.9 X 104 3.9 X 104 
10-20 2.9 5. 5.2 
20-30 3.1 4.5 4.5 
30-40 2:6 5.2 5.2 
40-50 2.2 5.2 5.2 
50-60 2.9 3.9 3.9 
60-70 3.1 2.4 2.4 
70-80 2.6 1.6 1.8 
80-90 7.8 — 2.4 
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Dilute Solution Viscosity of Polymers Prepared at 5°C. 

Another means of measuring molecular weight heterogeneity is by following 
the dilute solution viscosity of the polymer at various conversions. The 
viscosity changes with 75/25 and 66/34 charge ratios are shown in Fig. 3. It 
will be noted that with fresh latex the viscosity increases regularly with con- 
version but after the latex is aged the viscosity of the initial polymer has 
increased and that of the higher conversion samples has remained constant. 
The absence of a maximum as shown in Table I may be attributed to the larger 
amount of mercaptan present or to the lower polymerization temperature. 
Previous work (23) appears to have been done on aged latex. Similar data 
were collected for 50/50 and 25/75 butadiene to acrylonitrile charge ratio 
and are in Fig. 4. The viscosity increases to a maximum and then decreases 
for the 50/50 charge ratio. The decrease may be due to the increased rate 
of mercaptan disappearance at higher conversion together with the reduced 
tendency for branching of the substantially vinyl polymer. It seems that the 
dimethylformamide solvent consistently yields lower results than dichloro- 
benzene—isopropanol mixture. Again the viscosity of the aged latex has changed 


_ in the lower conversion ranges but has not increased above that of unaged 


latex in the higher conversion ranges. The few data with the 25/75 charge 
ratio indicate a regular decrease in viscosity which might be expected from 
the increasing acrylonitrile content of the polymer in addition to changes in 
the modifier disappearance curve and the reduced chances of branching. 





4 T by T i T T T .. * 
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Ye CONVERSION 
Fic. 3. Viscosity conversion curves. 75/25 butadiene/acrylonitrile charge ratio in benzene- 


isopropanol (e); 66/34 same (m); 66/34 in o-dichlorobenzene-isopropanol (a), same after two 
months (¥). 
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Fic. 4. Viscosity conversion curves. 50/50 butadiene/acrylonitrile charge ratio in dichloro- 


bensene-tsopropanol (e); same in dimethylformamide (m); same after two months (w); 25/75 
charge ratio in same (a). 


The increment viscosity may be calculated from the viscosity conversion 
data. The increment viscosity is that which will change the viscosity of the 
whole polymer from the value at the beginning of the range to the value at 
the end of the range. The data are in Table VI. The increment viscosity 
for the 75/25 and 66/34 ratios increases regularly, for the 50/50 charge ratio 
rapidly and then diminishes, and for the 25/75 ratio diminishes regularly. 


Viscosity Molecular Weights 


Conversion of the viscosity data into molecular weights can be made by 
using a Staudinger type equation. Almost any arbitrary equation would 


TABLE VI 
CALCULATED INCREMENT VISCOSITY FROM VISCOSITY—CONVERSION DATA 








75/25 charge ratio 66/34 50/50 25/75 
Conversion range 





Increment viscosity 








10-20% 0. 32 0.68 — a 
20-30 0. 68 0.95 1.83 0.55 
30-40 0. 68 1.28 3.31 0.15 
40-50 1.05 1.87 4.54 —ve 
50-60 1.57 1.92 2.50 —ve 
60-70 2.14 2.18 0.80 —ve 
70-80 a= 2.61 0 — 
80-90 = 2.82 0 — 
90-100 = 3.33 0 os 
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suffice for relative values. One applicable in this work has been published (18). 
The difference between the data for benzene and toluene solution is not 
sufficient to exclude the use of the equation for toluene solution when benzene 
or dichlorobenzene is the solvent. The intrinsic viscosities in all these solvents 
appear to be nearly identical in practice as shown in part of Fig. 3. The 
equation is 7 = 4.9 X 104.M*. This was used to calculate the molecular 
weight versus conversion, Table VII. The increment viscosity molecular 
weights were calculated likewise and are in Table VIII. 


TABLE VII 


CALCULATED VISCOSITY AVERAGE MOLECULAR WEIGHTS 























% Conversion 75/25 66/34 50/50 25/75 
charge ratio charge ratio charge ratio charge ratio 
10 <1x10° | <1 x 10 — — 
20 1.5 4.5 6.9 X 104 24.0 X 104 
30 3.0 | 8.0 18.5 18.5 
40 4.9 12.5 34.0 12.5 
50 7.8 | 17.8 45.0 6.9 
60 12.0 | 21.8 48.0 3.1 
70 17.5 24.5 43.0 <i 
80 — 29.5 38.0 — 
90 _— 34.5 34.5 —_ 
100 — | 37.5 29.5 _ 
TABLE VIII 


CALCULATED INCREMENT VISCOSITY MOLECULAR WEIGHTS 











Conversion 75/25 66/34 50/50 25/75 
range charge ratio charge ratio charge ratio charge ratio 
10-20 2.8 X 104 7.8 X 10! _— —_ 
20-30 7.8 18 3.9 X 10+ 7.2 X 104 
30-40 7.8 27 66 21 ; 
40-50 21 34 84 —ve 
50-60 34 40 53 . —ve 
60-70 46 46 14 —ve 
70-80 — 58 0 — 
80-90 — 58 0 — 
90-100 — 66 0 — 

















Heterogeneity of Polymers 
Having both the number average molecular weights and the viscosity aver- 
age molecular weights it is possible to estimate the heterogeneity of the 
polymer (8, 16). From the relationships and the viscosity molecular weight 
equation it can be derived (following the suggestion of Meehan*) that 
8 = 2.37 Vlog Rand My = M, R-*-*, 
where R is M./Mn. Using the above values, 8 and Mp) were calculated for 


* Private communication to participants in The Rubber Reserve research program. 
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the 66/34 and 50/50 charge ratios. The data are in Table IX. For both 
ratios the heterogeneity increases with conversion, and the peak of the distri- 
bution curve shifts to lower molecular weights. 

TABLE IX 


CALCULATION OF THE HETEROGENEITY AND PROBABLE MAXIMUM IN THE DISTRIBUTION 
CURVE OF MOLECULAR WEIGHTS 




















66/34 charge ratio | 50/50 charge ratio 
% Conversion 

B Mo | B | Mo 
10 —ve 94X10 | — — 
20 0 4.5 | 0.47 | 3.4.x 10 
30 0.59 3.2 1.46 1.9 
40 1.12 2.4 2.04 1.4 
50 1.35 2.2 2.30 | 1.2 
60 1.60 2.0 2.41 1.1 
70 1.84 1.4 2. 42 | 1.0 
80 | 2. 22 0.9 2.46 0.8 
90 | — | — 2.38 | 0.8 








In a similar manner the heterogeneity and the peak of the distribution curve 
for the increment polymer can be calculated, and the data are in Table X. 
The form of the curves is similar although the heterogeneity of the 50/50 ratio 
copolymer seems to be diminishing with increasing conversion. This again 
occurs when polymer of lower molecular weight is formed. 


From these calculated curves the usual derived data can be obtained. These 
include the percentage of polymer having various characteristics, etc. In 
general the values of 8 and Mp are sufficient to visualize the approximate 
molecular weight distribution of the polymers, and the data show that the 
system behaves much like GR-S or Buna S-3 (3) except as will be shown 
later when polymer rich in acrylonitrile is formed. 


Bound Acrylonitrile versus Conversion 
In addition to molecular weight heterogeneity there is chemical hetero- 
geneity along the chain. One aspect of this problem which can be studied is 


TABLE X 


CALCULATION OF HETEROGENEITY AND PROBABLE MAXIMUM IN THE DISTRIBUTION 
CURVE OF MOLECULAR WEIGHTS FOR THE INCREMENT POLYMER 

















Conversion 66/34 charge ratio 50/50 charge ratio 

range 

B Mo B Mo 
10-20 1.0 4.1 X 10! — — 
20-30 1.8 1.9 2.2 1.3 X 10! 
30-40 2.0 1.9 3.9 ‘3 
40-50 2.3 1.6 | 2.6 1.0 

2.4 0.9 2.4 1.0 
60-70 2.7 0.4 1.8 1 | 
70-80 2.9 0.2 | ae ai 

| 
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the variation of bound acrylonitrile with conversion. The -results for the 
75/25 butadiene/acrylonitrile charge ratio are in Fig. 5. The bound acrylo- 
nitrile values were determined experimentally and used to calculate the incre- 
ment acrylonitrile. This is the acrylonitrile content of the polymer being 
produced over a narrow range near the conversion indicated. Also the bound 
acrylonitrile was used to calculate the charge ratio of the residual monomers 
which would be expected to be related to the increment acrylonitrile through 
the reactivity ratios of the respective monomers. It is evident that the bound 
acrylonitrile is higher in the initial polymer than it is in the charge ratio and 
that polymer increasingly rich in butadiene is formed during the reaction to 
bring the final average polymer to a bound acrylonitrile equivalent to the 
charge ratio. During the reaction also the percentage of butadiene in the 
residual monomers increases. 
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Fic. 5. Acrylonitrile versus conversion for 75/25 butadiene/acrylonitrile charge ratio. (e) in 
polymer, (ms) in increment polymer, (a) in residual monomer. 


A different picture is observed when the 50/50 butadiene/acrylonitrile 
charge ratio is used, Fig. 6. These data show that the polymer formed is 
remarkably uniform chemically to 90% conversion after which there is a short 
period in which acrylonitrile content in the polymer formed rises steeply. The 
residual charge ratio increases slowly at first and then more rapidly until it 
appears as if all the butadiene is exhausted near the end of the reaction. 


Still a different picture is evident from the 25/75 butadiene/acrylonitrile 
charge ratio, Fig. 7. In this case the bound acrylonitrile starts at a value 
much lower than the charge ratio and increases slowly at first and then rapidly 
until the final average is equivalent to the charge ratio. The increment 
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Fic. 6. Acrylonitrile versus conversion for 50/50 butadiene/acrylonitrile charge ratio. (¢) in 
polymer, (wa) in increment polymer, (a) in residual monomer. 
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Fic. 7. Acrylonitrile versus conversion for 25/75 butadiene/acrylonitrile charge ratio. (e) in 
polymer, (¥) in increment polymer, (wa) in residual monomer. 


acrylonitrile likewise increases slowly at first and then rapidly. It attains 
the impossible values of slightly over 100%, which indicates that the slope of 
the bound acrylonitrile curve is changing too rapidly probably through an 
analytical error. Not only are bound acrylonitrile analyses difficult in this 
region but also the latex suddenly forms a suspengel which must be sampled 
by breaking the bottle and chipping off a portion of the spongy mass. This 
is an uncertain way of obtaining a representative sample. Also the mechanism 
of formation of this suspengel may be such that abnormal bound acrylonitrile 
values would result. This same error appears naturally in the residual charge 
ratio which indicates the disappearance of butadiene and its later reappearance. 
Both anomalies can be removed by slight changes in the original bound acrylo- 
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nitrile curves. Similar anomalies were observed in systems rich-in styrene (13) 
although the deviations were less pronounced and were definitely due to 
analytical errors in the region investigated. 


Reactivity Ratios of the Monomers, Q and e 


By extrapolating the bound acrylonitrile versus conversion curve to zero 
conversion an estimate of the reactivity ratio 7; and rz (where 7; is the ratio of 
the rate of reaction of butadieny] radical with butadiene to the rate of reaction 
of butadienyl radical with acrylonitrile and r2 is the ratio of the rate of reaction 
of acrylonitrylyl radical with acrylonitrile to the rate of reaction of acrylo- 
nitrylyl radical with butadiene) may be made by the linear method (2). The 
results are r; = 0.29 and r2 = 0.02. The value of 7; is lower than the previ- 
ously tabulated results (24) owing to the wider molar ratios used and possibly 
to the method of calculation which eliminates the undue weighting of the 
results obtained with small amounts of acrylonitrile in the charge. The value 
of r2 reported earlier (24) was —0.1, which was meaningless. The present 
value of r2 = 0.02 seems to fit the experimental results satisfactorily. 


It is interesting to plot the per cent acrylonitrile in the charge against the 
per cent acrylonitrile in the initial polymer formed. In this way an S-shaped 
curve is obtained, Fig. 8, in which the points are from Figs. 5, 6, and 7, and the 
curve is the theoretical for 7; = 0.29 and 72 = 0.02 in the equation 
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Fic. 8. Acrylonitrile in polymer versus acrylonitrile in charge. (e) data from Figs. 5,6, and 7. 
Curve theoretical for r, = 0.29 and rz = 0.02. 
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in which r is the molar ratio of butadiene to acrylonitrile in the monomers 


dB. ‘ : le ; 
and aa is the molar ratio of butadiene to acrylonitrile in the polymer being 


formed. For convenience these were calculated in terms of per cent acrylo- 
nitrile. The closeness of the agreement is satisfactory except for the charges 
very low in acrylonitrile. Were the copolymerization random, the result 
would be a straight line connecting zero and 100%. With butadiene and 
acrylonitrile, the polymer formed when less than 44% acrylonitrile is in the 
charge is richer in acrylonitrile than is the charge ratio, whereas if the acrylo- 
nitrile in the charge is greater than 44% the acrylonitrile in the polymer is 
less than in the charge ratio. This illustrates the strong tendency for butadiene 
and acrylonitrile to alternate in the chain, an effect which will be discussed 
later in more detail. 


The data of Smith (21) may be used to calculate the distribution of acrylo- 
nitrile between an oil phase and water, assuming that the oil phase behaves 
similarly to styrene. This may not be true but until a satisfactory study* of 
the partition of acrylonitrile between butadiene and water is available the 
data (21) are at least an approximation. The calculated values used in the 
study are in Table XI. 


TABLE XI 


CALCULATION PERCENTAGE OF THE ACRYLONITRILE IN THE OIL PHASE 











Butadiene /acrylonitrile charge ratio Acrylonitrile % in oil phase 
90/10 56 
80/20 62 
70/30 65 
60/40 69 
50/50 71 
40/60 72 
30/70 73 
20/80 74 
10/90 75 
0/100 76 








Recalculating the results in Figs. 5, 6, and 7 on this basis yields values of 

= 0.18 and r. = 0.03. Plotting the data in Fig. 9 shows that now the 
experimental results are close to the theoretical curve for r; = 0.18 and rz = 
0.03 all along the line. At least a satisfactory approximation has been achieved. 
* A preliminary study of the distribution of acrylonitrile between butadiene and water suggests 

that under the experimental conditions used in this paper 87%, 92%, and 94% of the acrylonitrile 
is in the oil phase for 75/25, 50/50, and 25/75 butadiene /acrylonitrile charge ratio respectively. 


These results applied to Figs. 5,6, and 7 yield r, = 0.24 and r2 = 0.03. Such data would yield 
points and a theoretical curve similar to but intermediate between those of Figs. 8 and 9. 
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Fic. 9. Acrylonitrile in polymer versus acrylonitrile in oil phase of charge. (e) data from 
Figs. 5, 6, and 7 as changed by Table XI. Curve theoretical for r, = 0.18 and r. = 0.03. 





Discussion 


A general comparison of interest is the increment acrylonitrile curves. A 
series has been published by Semon (19, 20). His results showed that 62.6/ 
37.4 mole ratio of butadiene to acrylonitrile was azeotropic, i.e., the polymer 
was of constant composition throughout. Also his data showed that when 
the charge ratio varied from 56/44 to 77/23 the initial polymer formed varied 
from 1.5/1 to 2/1 mole ratio of butadiene to acrylonitrile. The experimental 
data above are qualitatively similar but suggest an azeotropic mixture con- 
taining a higher acrylonitrile proportion. The exact value would seem to be a 
fortuitous balance between reactivity ratio, distribution of monomers between 
phases, and perhaps other factors, and it would be surprising if these could be 
balanced exactly to 100% conversion. Also the data above suggest that the 
initial polymer produced at 5°C. is richer in acrylonitrile than for the corre- 
sponding charge ratio at higher temperatures since the mole ratio of the initial 
polymer when the charge ratio is varied from 25/75 to 75/25 varies from 1/1 
to 1.75/1 mole ratio butadiene to acrylonitrile, i.e., lower than at the higher 
polymerization temperature. This may be due to lower solubility of acrylo- 
nitrile in water at the lower temperatures. The azeotropic charge ratio from 
Fig. 9 would appear to be 54/46 butadiene/acrylonitrile for bulk phase 
polymerization. Comparison with similar curves for the styrene—acrylonitrile 
system (21) shows the more flattened central portion of the curve and the shift 
upwards of the azeotropic concentration of the above data. 


The dilute solution viscosity indicates that the molecular weight increases 
usually regularly but in the case of the 50/50 ratio to a maximum and then 
decreases and in the case of the 25/75 ratio polymer decreases regularly. In 
the former the changed direction of the trend coincides with the formation of 
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low molecular weight polymer owing to increased rate of disappearance of the 
mercaptan, and in the latter to this effect and the formation of polymer richer 
than 50% in bound acrylonitrile. The high acrylonitrile polymer seems to 
behave in a different fashion from the less rich copolymer in the same solvent, 
perhaps in a manner similar to a polyelectrolyte but under the influence of 
those forces associated with the oil resistance of this type of polymer. 


The relationship between viscosity and molecular weight is sufficiently close 
to be of suitable use for calculation of approximate viscometric molecular 
weights. These compared with the number average molecular weights yield 
some indication of the probable maximum of the molecular weight distribution 
and the heterogeneity of the various polymers. Since a number of effects 
have to be considered negligible the results can be in considerable percentage 
error but probably little relative error. In general as the conversion increases 
the heterogeneity becomes greater and the probable maximum in the distri- 
bution curve occurs at lower molecular weights. However, when the viscosity 
curve reverses direction with the formation of low molecular weight polymer 
the relationships again do not follow the expected pattern, i.e., the hetero- 
geneity passes through a maximum and the peak of the distribution curve 
remains relatively constant. The data obtained when polymer containing 
over 50% acrylonitrile is formed are difficult to interpret as yet. 


It will be noted that in copolymerizations with butadiene or styrene with 
acrylonitrile (11) there is a greater tendency for these monomers to form a 
radical more reactive towards acrylonitrile than to themselves, whereas acrylo- 
nitrile also forms a radical and more reactive towards butadiene or styrene 
than to acrylonitrile. Such a situation would cause the polymer chain to 
tend to alternate between acrylonitrile and butadiene or styrene. Thus the 
rate of polymerization is greatest for approximately equal molar ratios. This 
also leads to a small value of the 7,.72 product (10) which is considered a 
measure of the alternating tendency. Price (17) has used such data also to 
compute a value Q and a value e, the former to represent the average reactivity 
of the double bond and the latter to represent a measure of the polarity of the 
complex. The Q and e for acrylonitrile in the butadiene—acrylonitrile system 
appear to be 0.74 and 1.47 respectively. The values are appreciably higher 
than those recorded (see (11)) owing to the greater tendency of butadiene and 
acrylonitrile to alternate under the conditions of these experiments. This has 
been observed for the butadiene-styrene system (16). Further investigation 
under more varied conditions is needed to ascertain whether such effects are 
real. 
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THE COMPOSITIONAL HETEROGENEITY OF 
BUTADIENE-STYRENE COPOLYMERS 
SYNTHESIZED AT —18°C. IN EMULSION! 


By R. J. ORR AND H. LEVERNE WILLIAMS 


Abstract 


A study has been made of butadiene-styrene copolymers formed at —18°C. 
From analyses for bound styrene in the product for various conversions and 
initial butadiene-styrene ratios the reactivity ratios were calculated to be 
r, = 1.37 and rz = 0.38 compared with 1.8 and 0.6 at 45°C. Q and e for 
butadiene were 1.38 and 0.008 relative to styrene at 1 and —0.8. Increment 
bound styrene curves calculated for each stage of the reaction indicated that the 
polymers were remarkably homogeneous at low conversions. The chain transfer 
reaction using mixed tertiary mercaptans as the modifier was studied. Regu- 
lating indices were found to have decreased with polymerization temperature. 
Number average (./,,) and viscosity average (My) molecular weights were calcu- 
lated from mercaptan disappearance and vistex intrinsic viscosity data re- 
spectively. The molecular weight heterogeneity increased with increasing con- 
version and initial mercaptan content. The increment number average molecular 
weights were found to diminish with conversions, whereas the increment viscosity 
average increased at higher conversions as conversion increased. 


Introduction 


While general purpose chemical rubbers are products containing about 23% 
bound styrene it is possible to copolymerize butadiene and styrene in all ratios 
to copolymers which show a continuous variation in their properties between 
those of polybutadiene and polystyrene. Recently it has been shown that 
polymerization temperatures much lower than formerly used yielded a superior 
product when butadiene was present as a major component of the monomer 
blend (31, 32). Among reasons for this improvement have been advanced 
(3, 4, 8, 13, 14, 20, 25, 33) the suggestions of greater regularity of the chain and 
the formation of a more linear polymer chain structure. 


The preparation of copolymers using the present techniques gives rise to 
products showing chemical and physical heterogeneity. The chemical hetero- 
geneity results from variations in the proportions of the two copolymerized 
species along the polymer chain. This is a result of the difference in activity 
of the two monomers which governs the relative rates with which the two 
monomers enter the copolymer. The physical heterogeneity is a result of 
the presence of a range of molecular weights in the product. The distribution 
of the molecular weight is for the most part a function of the types of chain 
transfer agents present and their mode of disappearance. While all these 
factors are well known for the butadiene-styrene copolymerizations at 50°C.., 
it has become most desirable to know how they change with temperature. 


1 Manuscript received in original form September 26, 1950, and, as revised, November 28, 
1950 
Contribution from the Research and Development Division, Polymer Corporation Limited, 
Sarnia, Ontario. Presented in part before the Chemical Engineering Subject Division, Chemical 
Institute of Canada, Toronto, June, 1950. 
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Techniques 


The polymerizations were carried out in 8 oz. peroxide bottles with punctured 
metal screw caps (10) in which the fiber gasket was replaced by a double rubber 
gasket consisting of an outer oil resistant Koroseal* sheet and an inner self 
sealing Butyl** gasket prepared from undercured unpigmented sheets. 


The recipe used in the synthesis of the polymer is shown in Table I. The 
development of this recipe was reviewed earlier and the techniques outlined in 
detail (28). 








TABLE I 

CHARGE RECIPE FOR —18°C. POLYMERIZATION 
Monomers 100 
Potassium fatty acid S.F. soap flakes (Swift Co.) 5.1 pt./100 monomers 
Daxad 11*** 0.1 
Methanol 46 
Water, distilled 157 
MTM**** As stated 
FeSO,. 7H:O0 0.20 
As.O3 0. 034 
Sequestrene AA***** 0.12 
K,4P207 0. 056 
Diisopropylbenzene monohydroperoxide (Hercules) 0.15 
Ditertiary butyl hydroquinone 0.2 








* A product of B.F. Goodrich Chemical Co. 
** A product of Polymer Corporation Lid. 
*** Trade name of a dispersing agent, Dewey and Almy Chemical Co. 
*#4* 4 mixed tertiary mercaptan consisting of 60% Ci2, 20% Cis, 20% Cis, average molecu- 
lar weight 220, Phillips Petroleum Co. 
er Ethylene dinitrilo tetraacetic acid, Alrose Chemical Co. 


After the reaction was stopped the degree of conversion of monomer to 
polymer was determined (12) and the unreacted monomers were removed from 
the latex by low pressure steam. The latex was filtered through glass wool to 
remove precoagulum and coagulated in isopropanol. The polymer was washed 
with water, dried, and milled or pressed to a thin sheet suitable for determina- 
tion of refractive indices. The refractive index of this sheet was determined 
at 60°C. (26) and the styrene content of the polymer calculated. Nocorrection 
was made for density changes of the polymer due to its more regular structure 
when produced at lower temperatures (25). The mercaptan in the latex was 
determined by amperometric titration (15, 16). The intrinsic viscosities of 
the polymers formed were measured at 29.9°C. by dissolving the latex directly 
in a 80:20 benzene-isopropanol solution. Dilutions were made with pure 
benzene so that intrinsic viscosities in benzene could then be calculated (7, 
11, 26). 

Experimental 
Bound Styrene versus Conversion 


The per cent styrene in the polymer was measured for various initial ratios 
of butadiene to styrene in the charge and for various conversions. The data 
are in Fig. 1. In general there is an appreciable conversion range over which 
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Fic. 1. Bound styrene versus conversion. The charge ratios top to bottom are 30/70, 50/50, 
60/40, 70/30, 72/28, and 75/25 butadiene to styrene. 


the bound styrene does not change much before it rises to correspond to the 
charge ratio. The length of this plateau also tends to decrease with increasing 
styrene content in the charge. Somewhat similar results have been reported 
elsewhere for single charge ratios at 5°C. and —20°C. (25, 33). Additional 
data obtained in a comparable study at a single charge ratio (X—435 recipe) 
at 5°C. are in Table II. Comparison of these values with those obtained at 


TABLE II 


BOUND STYRENE VERSUS CONVERSION AT 5°C. 
FOR 72/28 CHARGE RATIO (27) 











Conversion, % % Styrene in polymer 

5.6 20.7 
14.1 20.0 
24.3 20.1 
28.8 21.1 
34.3 20.5 
39.8 20.7 
49.6 21.7 
49.8 21.3 
52.6 21.7 
62.6 21.2 








higher temperatures (21, 22, 24, 25) shows that the tendency towards constant 
bound styrene at low conversions is characteristic of low temperature polymer- 
izations. However, there is a comparatively small change in the chemical 
heterogeneity with temperature. 


Extrapolation of the curves in Figure 1 to zero conversion yields the data in 
Table ITI. 
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TABLE III 


BOUND STYRENE AT ZERO CONVERSION BY EXTRAPOLATION 
OF THE DATA IN FIG. 1 











Charge ratio % bound styrene 
30/70 57.0 
50/50 41.0 
60/40 30.8 
70/30 22.8 
72/28 21.2 
75/25 18.1 








Monomer Reactivity Ratios, Q and e 


Use of the linear method for calculating reactivity ratios (6) yielded 7; = 1.37 
and r2 = 0.38, where 7; is the ratio of the rate of reaction of butadienyl radical 
with butadiene to the rate of reaction of butadienyl radical with styrene, and 
r. is the ratio of the rate of reaction of styryl radical with styrene to the rate 
of reaction of styryl radical with butadiene. Alfrey and Price (2) describe 
- the behavior of monomers in terms of two factors, Q the average reactivity 
and e, the polar properties of the monomers. For this treatment styrene was 
taken as the reference monomer having a Q of 1 and an e of —0.8. If the 
reactivity ratios did not change appreciably with temperature, butadiene 
would possess the same Q and e values in relation to styrene as it did at 50°C., 
i.e., 1.33 and —0.8 (29). The values obtained at —18°C. are 1.38 and 0.008 
respectively. While Q has remained unchanged e has changed greatly. 
€ appears very sensitive to changes in 7; or r2. In the consideration of a large 
number of monomer reactivity ratios, Lewis, Mayo, and Hulse (18) pointed 
out that r; X re equals unity if the monomers are entering the polymer chain 
completely at random. In the present system 7 X r2 is 0.52 which indicates 
a slight tendency to alternation in the butadiene-styrene copolymer. 


Increment Bound Styrene, Residual Unreacted Styrene 


The increment bound styrene is the average composition of the polymer 
formed over an increment in conversion. In Table IV are the increment bound 
styrenes for the various charge ratios. It may be seen that for 30/70 charge 
ratio the increments are quite constant to 30-40% conversion and for the 
remaining charge ratios the increments are quite constant up to about 50-60% 
conversion. After this period the polymer increases in styrene content as the 
conversion proceeds usually rapidly at first and then more slowly. It is evident 
from Table II that similar results were obtained at 5°C. (27). 


The residual styrene in the monomers may be calculated at various con- 
versions. The data are in Table V. The residual styrene in the charge rises 
with conversion in a manner very similar to the bound styrene of the increment 
polymer (Table IV). The two values should be related through the reactivity 
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TABLE IV 
INCREMENT BOUND STYRENE CALCULATED FROM FIG. 1 




































































% bound styrene in increment polymer for 
Conversion 
range, % 30/70 | 50/50 60/40 | 70/30 | 72/28 | 75/25 
charge | charge charge | charge charge charge 
ratio ratio ratio ratio ratio ratio 
0-10 57.6 41.6 31.0 22.9 21.3 18.6 
10-20 59.0 42.8 31.8 23.1 21.7 19.8 
20-30 60. 4 44.0 32.3 23.6 21.8 20.7 
30-40 64.6 45.2 32.9 23.6 22.0 21.7 
40-50 73.6 46. 4 33.0 23.8 22.2 22.7 
50-60 75.8 57.2 32.6 24.6 27.0 23.1 
60-70 73.0 60. 2 45.2 32.0 23.4 | 24.6 
70-80 76.0 59.4 52.8 37.6 25.4 | 26. 4 
80-90 78.4 53.2 57.4 | 39.9 37.0 33.0 
90-100 80.8 | 50.0 49.0 | 48.9 61.8 | 39.0 
TABLE V 
RESIDUAL STYRENE IN CHARGE VERSUS CONVERSION 
% Residual styrene in charge for 
Conversion, %| 30/70 50/50 60/40 70/30 72/28 | 75/25 
charge charge charge charge charge | charge 
ratio ratio ratio ratio ratio | ratio 
0 70 50 40 30 28 25 
10 71 51 41 31 28 26 
73 52 42 32 36 26 
30 75 53 43 33 31 27 
40 76 54 45 34 32 29 
50 77 56 48 37 34 29 
77 56 51 39 37 31 
70 78 54 53 42 42 | 33 
79 52 53 44 49 | 36 
90 81 50 49 49 62 39 
{ 
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Fic. 2. Per cent styrene in polymer versus per cent styrene in monomer blend. (e) calculated 
from experimental data. The curve is for r, = 1.37 and re = 0. 38. 
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ratios r; and rs. When the residual styrene in the charge is plotted against 
the styrene content of the increment polymer the points in Fig. 2 are obtained. 


The line is the theoretical for 7, = 1.37 and r2 = 0.38 using the relationship 


B 
dB = B[ n\S) +1 


dS S$ B 
a + ro 


where S is the concentration of styrene and B is the concentration of butadiene. 
In order to make the information of more practical use the calculations were 
done using the per cent styrene in the increment polymer and in the residual 
monomers rather than mole ratios. Although the agreement may at first 
sight appear poor, the results are actually good. The values of 7; and rz were 
calculated from low conversion data, and the points obtained at low conversion 
lie on the theoretical curve. Points obtained at high conversion lie below the 
line for high butadiene charges and above the line for high styrene charges, 
that is, 7; and 72 both increase with conversion. The data attained with the 
50/50 charge ratio were off the curves both as used to calculate 7; and 72 and in 


. Fig. 2. The cause for this anomaly is not known. 


Modifier Disappearance Curves 


Another form of heterogeneity is that of molecular weight. The method of 
investigation of this phase of the problem was by a study of the chain transfer 
mechanism at this temperature together with the measurement of viscosity 
average molecular weights. Three different initial mercaptan charges were 
used with a 70/30 charge ratio and two other charge ratios were investigated 
similarly with one modifier content. The data on mercaptan disappearance 
are in Fig. 3. It is evident that the mercaptan disappears slowly and is not 
markedly affected by concentration of mercaptan or of monomer charge ratio. 








LOG % RESIDUAL MERCAPTAN 








% CONVERSION 


Fic. 3. Modifier disappearance curves. The charge ratios butadiene to styrene and parts 
mercaptan are upper (top to bottom) 70/30 and 0. 12, 70/30 and 0. 24, 70/30 and 0. 48; and lower 
(top to bottom) 50/50 and O. 24 (ms), 70/30 and 0. 24 (a), and 30/70 and 0. 24 (e). 
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If the natural logarithm of the residual mercaptan is plotted against the 
per cent conversion, it is found that a linear relationship exists at low con- 
versions. The slope of the linear portion has been adopted as a measure of 
the relative rates of the chain transfer reaction and conversion of monomer 
to polymer. These slopes were determined from Fig. 3 and the results are 
given in Table VI. The accuracy of the method is not great so that the results 
are not considered exact. 


TABLE VI 


REGULATING INDICES OF MTM at —18°C. 











Charge ratio Initial RSH content Regulating index 
parts/100 monomers 
70/30 butadiene/styrene 0.12 0.5 
70/30 = 0.24 0.7 
70/30 si 0.48 1.4 
50/50 = 0. 24 0.5 
30/70 ii 0.24 0.6 
72/28 e (27) 0. 20-0. 60 7.8 











It may be seen from this that the regulating index under all conditions is 
quite small. There appears to be a trend toward larger regulating indices at 
higher initial mercaptan concentrations. Comparison with higher tempera- 
ture (7, 26) fatty acid soap systems where values in the range of 1.5 to 2.5 
were obtained (compared with 7.8 at 5°C. (27) using a rosin acid soap) indi- 
cates that the chain transfer reaction slows down more with temperature than 
does the rate of propagation. The effect of the methanol antifreeze used at 
— 18°C. on the modifier disappearance curves has yet to be determined. 


Number Average Molecular Weights 

Assuming that for each molecule of mercaptan consumed there was one 
polymer molecule formed, number average molecular weights can be calculated 
by the formula: 


Number average molecular weight (M,) = Total weight of polymer. 
Number of moles of polymer 





The results are in Table VII. The assumption is valid at least until branching 
or cross linking begins (see below). Termination reactions by other than 
mercaptans tend to decrease the number average molecular weights. Thus 
the values calculated for low conversions are subject to minimal error. 


The results show that the number average molecular weight decreases with 
increasing modifier in the charge and with increasing conversion. The latter 
effect is due to the increased rate of consumption of mercaptan at the higher 
conversion, Fig. 3. At 5°C. (27) the reverse is true initially since the rate of 
modifier disappearance diminishes with increasing conversion then increases; 
the number average molecular weight increases with conversion. These 
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TABLE VII 


CALCULATED NUMBER AVERAGE MOLECULAR WEIGHTS 














50/50 30/70 72/28 
70/30 charge ratio charge charge _ |charge ratio 
Conversion, % ratio, ratio, at LU. 
0.12 part | 0.24 part | 0.48 part | 0.24 part | 0.24 part | 0.40 part 
MTM MTM MTM MTM MTM MTM (27) 
10 37 X 10‘ | 13 K 104 | 3.5 KX 104} 18 K 104 | 18 X 104 | 1.1 X 104 
20 37 14 4.0 18 15 i 
30 39 14 3.9 18 13 2.2 
40 31 13 3.8 15 11 2.8 
50 25 9.5 3.9 13 9.8 3.6 
60 22 8.1 4.1 13 8.9 6.1 
70 20 42 4.2 11 8.2 — 
80 — — 4.4 11 —_ — 
90 —_ _— 4.6 ll —_ — 























results are illustrated more clearly by the increment number average molecular 
weight. These are the molecular weights of the polymers formed over the 
conversion range indicated. The data are in Table VIII. 


Intrinsic Viscosity versus Conversion 


The intrinsic viscosity versus conversion curves are in Fig. 4. Inspection 
of the intrinsic viscosity data indicates that the viscosity—conversion curve is 
quite flat until relatively high conversions are reached. No gel-point was 
apparent. This indicates that at these lower temperatures gel formation and 
hence the cross-linking reaction is greatly suppressed, leading to formation of 
a greater predominance of long linear polymers. The intrinsic viscosity did 
show a rise at higher conversions. Both this rise and the plateau preceding it 
were functions of the initial mercaptan content in the manner shown in 


Table IX. 


TABLE VIII 


CALCULATED INCREMENT NUMBER AVERAGE MOLECULAR WEIGHTS 














70/30 charge ratio 50/50 30/70 72/28 
Conversion charge ratio,|charge ratio,|charge ratio 
range, % 0.24 part | 0.24 part jat 5°C. (27), 
0.12 part | 0.24 part | 0.48 part MTM MTM 0. 40 part 
MTM MTM MTM MTM 
0-10 3.7 X 104 | 1.3 X 104 | 0.35 X 104| 1.8 XK 104 | 1.8 X 104 | 1.1 x 104 
10-20 3.7 1.5 0. 46 1.8 1.3 3.0 
20-30 3.7 1.5 0.38 1.8 0.92 6.5 
30-40 2.0 0.92 0.35 0.92 0.83 12.2 
40-50 1.5 0.48 0. 42 0.83 0.65 11.0 
50-60 1.3 0.46 0.51 0.83 0.61 9.2 
60-70 1.2 0. 54 0. 57 1.0 0.59 - 
70-80 — — 0.61 0.92 _ — 
80-90 _ — 0. 64 0.97 -- — 
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Fic. 4. Viscosity conversion curves. The charge ratios butadiene to styrene and parts mer- 
captan (top to bottom) are 70/30 and 0. 12, 70/30 and 0. 24, 50/50 and 0.24 (@), and 70/30 and 
0.48 (m), and 30/70 and 0. 24. 


TABLE IX 


INFLUENCE OF INITIAL MERCAPTAN CONTENT ON CHARACTERISTICS OF INTRINSIC 
VISCOSITY CONVERSION CURVES FOR 70/30 CHARGE RATIO 








Initial mercaptan content | Conv. % before showing vistex rise | Height of viscosity plat. 





0. 12 parts (0.00054 mole) | 45 | [n] = 4.5 
0.24 parts (0.00108 mole) | 80 | 3.1 
0.48 parts (0.00216 mole) | 90 2.2-2.3 





As the styrene content of the polymer is raised, the vistex viscosity decreases 
perhaps owing to the role of styrene as a chain transfer agent or to a change in 
the viscosity—molecular weight relationship. 


The viscosity conversion curves appear to be more similar to those obtained 
in the high styrene range at higher polymerization temperatures than to those 
obtained in the low styrene range (26). Somewhat intermediate results are 
obtained at 5°C. For example (27) the conversion and corresponding intrinsic 
vistex viscosity for 72-28 charge ratio and 0.40 part MTM modifier are 25.2%, 
0.49; 35.8%, 1.06; 47.8%, 1.18; 60.0%, 1.47; and 67.0%, 2.05. At 5°C., 
variation in the amount of modifier changes the form of viscosity—conversion 
curves in the usual manner, i.e., more modifier lowers and flattens the curve, 
less raises and humps the curve. At — 18°C. this effect, if present, is minimized. 


Viscosity Molecular Weight 


It is possible from intrinsic viscosity data to calculate viscosity—average 
molecular weights by means of the well known equation 7 = KM:*, where 7 
is the intrinsic viscosity, 
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M is the viscosity average molecular weight, and 
K and a are constants characteristic of the polymer, solvent, and tempera- 
ture at which 7 was measured. 


Almost any reasonable arbitrary curve would suffice for relative values but 
it was of interest to estimate a most probable one. At 65% conversion, K and 
a of polybutadiene in toluene at 30°C. show the following dependence on poly- 
merization temperature (13). 


Polymerization temperature K a 
50 (°C.) 72.5 X 104 0.45 
5 26.4 X 104 0.55 
—20 10.6 X 104 0. 63 


For polystyrene polymerized at 60°C. K is 1.6 and a is 0.69 in toluene at 
30°C. (1). This was assumed not to change for polymer prepared at still 
lower polymerization temperatures. For GR-S in toluene at 30°C. K is 5.25 
and a is 0.667 (30). Inspection of the data on butadiene indicates a nearly 
linear relation between a and 7, whereas K and T are only approximately 
linear. Assuming linearity then the plot K or a versus temperature for poly- 
butadiene and polystyrene intersect or if they do not in fact intersect they 
approach very close for polymers prepared at very low polymerization tempera- 
tures. At that polymerization temperature at which K and a for polybutadiene 
are equal to or approach K and a for polystyrene, all copolymers of these 
monomers.should have a K and a equal to or approaching that of polybutadiene 
and polystyrene. This polymerization temperature is of the order of —28°C. 
at which K = 1.6 and a = 0.69 in toluene at 30°C. for polybutadiene, poly- 
styrene and copolymers. Drawing a line from K and for GR-S to this inter- 
section it is possible to estimate K and a for polymer of the GR-S type prepared 
at —18°C. The equation is 7 = 2.0 X 10* M°® at 30°C. in toluene. Some 
calculated data are in Table X using this equation. It is also evident from the 
data above that at a polymerization temperature of 50°C., introduction of a 


TABLE X 


CALCULATED VISCOSITY AVERAGE MOLECULAR WEIGHTS 








50/50 30/70 72/28 














70/30 charge ratio charge charge charge ratio 
Conversion, % ratio, ratio, at 5°C. (27), 
0.12 part | 0.24 part | 0.48 part | 0.24 part | 0.24 part | 0.40 part 
MTM MTM MTM MTM MTM MTM*, 
10 21 X 10° | 12 X 10° | 7.7 XK 10° | 7.8 X 10° | 3.8 < 10° | 0.62 x 105 
20 21 12 7.9 7.8 3.8 1.08 
30 21 12 rae 7.8 3.8 1.30 
40 21 12 rs 7.8 3.8 1.30 
24 12 ree i 7.8 3.8 1.50 
60 29 13 at 7.9 3.8 1.98 
70 — 13 eS 8.2 3.8 | 4.0 
80 — 14 7.8 9.0 3.8 | — 
90 — 15 8.2 10.0 — — 

















* Assuming K = 3.0 X 10°. 
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relatively large amount of butadiene into the polymer causes only a very small 
changein Kanda. At —20°C. the separation between K and a for polystyrene 
and polybutadiene is much less. Thus it seems evident that K and a derived 
above should show little dependence on the bound styrene content of the 
polymer. In general the value cannot be sufficiently in error to invalidate the 
results for the purpose of this paper. 


Similarly the increment viscosity average molecular weight may be calcu- 
lated. The data are in Table XI. 


TABLE XI 


CALCULATED INCREMENT VISCOSITY AVERAGE MOLECULAR WEIGHTS 




















| | 
| 50/50 | 30/70 | 72/28 
Conversion | 70/30 charge ratio charge | charge |charge ratio 
range, % ratio, | ratio, jatec. G7), 
| 0.12 part | 0.24 part | 0.48 part | 0.24 part | 0.24 part | 0.40 part 
MTM MTM | MTM MTM | MTM | MTM* 
0-10 © |X 105 12 X 10° | 7.7_X 105 | 7.8 X 10° | 3.8 X 10° | 0.62 x 105 
10-20 2.1 12 oe . ee) ee ae * 
20-30 | 2.1 me ae 7.8 3.8 | 1.75 
30-40 | 2.1 ae ft oe 7.8 3.8 | 1.30 
40-50 3.8 oe. | ae 4 7.8 38 | 32.5 
50-60 5.7 S + Be} 8.4 3.8 | 4.8 
60-70 a: a 8.80 | 10.3 3.8 | 2.3 
70-80 — 21 Tet 14.5 a3 | — 
80-90 — 24 «| «(18.5 17.0 —- | = 





* Assuming K = 3.0 X 104. 


Discussion 


Of-considerable interest is the change in the copolymerization constants 
with temperature. The values of r; equal to 1.37 and r2 equal to 0.38 are 
quite definitely lower than the previous values for butadiene and styrene in 
emulsion measured at 45°C. (24, 26), i.e., 71 equal to 1.8 and re equal to 0.6. 
The decrease is in the same direction as other investigators have found (9, 19, 
23) and may mean that the energy of activation of homo-propagation is slightly 
more than for heteropolymerization. A similar conclusion was reached in a 
concurrent study (5) of the butadiene—acrylonitrile system. Calculation of 
r, and r2 by the method of Mayo and Lewis (cf. 23) yielded the values r; equal to 
1.66 and rz equal to 0.42. Since these values were obtained in the range of 
10-25% conversion the somewhat higher values are in the right direction but 
the increase is rather surprising. The results do not change the conclusions 
reached in this paper since the product of 7; and r2 is changed only to 0.70. 
Calculation of the Q and e values corresponding to r; = 1.37 and r2 = 0.38 
yields the interesting data that Q remains unchanged with polymerization 
temperature but e becomes markedly less negative. This is relative to styrene 
with the same values as are used at higher temperatures. Thus the question 
arises as to a suitable standard for e at various temperatures. It is interesting 
to note that the change in e appears greater than that predicted (24), i.e., that e 
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should vary as the reciprocal of the square root of the absolute temperature. 
A similar conclusion was reached in a concurrent study (5) of the butadiene- 
acrylonitrile system. This may mean that the radius of the intermediate 
complex has changed also with temperature since it is not considered likely 
that the dielectric constant has altered substantially under these conditions. 
The effect of any methanol in the oil phase would appear to be negligible 
also (19). 


An approximation to the molecular weight distribution (17) may be de- 
termined from the viscosity average molecular weight (M,) and the number 
average molecular weight (M,) by the logarithmic number distribution Max- 
wellian type with unlimited range toward infinite molecular weight, i.e., 

dn = Ne” dy 
Vr 
where y = 1/8 In M 
Mo ; 
dn = numer of molecules in the interval y to y + dy, 
N = total number of molecules, 
8 = degree of heterogeneity and is zero for a homogeneous polymer, and 
Mp, = molecular weight at which peak occurs in the distribution curve. 


From the above relationships it may derived for butadiene-styrene co- 
polymers (following the suggestion of E. J. Meehan) that: 


B = 2.35-Vlog R, and My = M,, R®-*8, where R = M,. 
M, 
Also, if y = 1/ In M/ Mb, it may be shown that F,,, the weight fraction of the 


polymer having less than the molecular weight M is equal to 





z 


1 
we | e~* dz, where z = y — 4. 


-© 

This is the familiar probability integral for which values have been tabulated. 
When F,, is plotted against M it gives the integral molecular weight distri- 
bution of the polymer. Some data are in Table XII. 


Comparison with 8 determined at higher temperatures (7, 6th polymer 
Table XII) reveals that GR-S at 20% conversion had a 6 of 1.30 when 0.5 
part of dodecyl mercaptan is used as regulator. With 0.48 part (0.00216 
mole) MTM and at 20% conversion 8 had a value of 2.65. This larger value 
of 6 may be due to decreasing polymerization temperature. This difference 
is much greater than any error introduced by the estimated viscosity molecular 
weight relationship. 


The values of 8 are probably maximal since errors in the choice of K tend to 
decrease 8. The large values of 8 reflect the very large ratio of My to Mn as 
calculated from the experimental data and may be a measure of the 
structural changes in the polymer formed at — 18°C. relative to that formed at 
50°C. It may also be related to the statistical assumptions which may no 
longer be valid when the chain transfer reaction is reduced. 
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TABLE XII 





RATIO OF VISCOSITY TO NUMBER AVERAGE MOLECULAR WEIGHT, DEGREE OF HETEROGENEITY, 
AND MOLECULAR WEIGHT AT WHICH THE PEAK IN THE DISTRIBUTION CURVE OCCURS 
























































70/30 charge ratio 70/30 charge ratio 70/30 charge ratio 
Conv., % 0.12 part MTM 0.24 part MTM 0.48 part MTM 
R | B | Mo R B Mo R | B | Mo 
| 
10 5.68| 2.04 | 1.3 X 10| 9.25| 2.31 | 0.34 x 105| 22.0 | 2.73 | 5.5 x 108 
20 5.68| 2.04 1.3 8.60/ 2.10} 0.38 | 19.24| 2.67 6.8 
30 5.48| 2.02 | 1.4 8.60) 2.10 0.38 |19.74/ 2.68) 6.5 
40 6.77| 2.14| 0.98 9.25| 2.31 0.34 | 20.5 | 2.69} 6.2 
50 9.61| 2.33| 0.64 | 12.65) 2.35 0.21 19.74| 2.68) 6.4 
60 | 13.18/ 2.58| 0.46 | 16.05) 2.58 0.15 | 18.76) 2.65| 7.1 
70 | | 17.35| 2.62| 0.14 |18.34| 2.64| 7.3 
80 17.72| 2.62} 7.9 
90 | | 17.83| 2.63 | 8.2 
| | 
| 72/28 charge ratio 
50/50 charge ratio, 30/70 charge ratio, at 5°C. (27), 
0.24 part MTM 0.24 part MTM 0.40 part MTM 
10 4.33| 1.87 | 7.5 X 10'| 2.12] 1.34 | 11.5 x 10*| 5.62] 2.04 | 3.9 x 108 
2 | 4.33] 1.87| 7.5 2.53/1.49| 86 | 6.34|2:10| 5.6 
30 4.33) 1.87 | 7.5 2.92| 1.60! 6.8 5.65| 2.04 | 3:3 
40 5.20) 1.99 | 5.6 3.45] 1.72| 5.3 | 7.22) 2.00; 8.5 
50 6.00| 2.07 | 5.0 3.87/ 1.80} 4.4 | 4.05) 1.83] 16.0 
60 6.08/ 2.08 | 4.4 4.27| 1.87 | 3.7 3.25| 1.68 | 30.0 
70 7.36/2.19| 3.3 4.65|1.92| 3.3 | 
80 8.08| 2.24; 3.2 | | 
90 | 9.00) 2.29 | 3.0 | | | | | 
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THE CALCULATION OF THE CHANGE IN VIBRATIONAL 
FREQUENCY DUE TO CHANGES IN MASS, GEOMETRY, 
AND POTENTIAL CONSTANTS OF A MOLECULE! 


By H. J. BERNSTEIN 


Abstract 


An expression for the change in vibrational frequency in terms of changes in 
the kinetic and potential energy of a molecule has been derived. 

For numerical calculations approximations are given. The change in fre- 
quency due to configurational changes has also been estimated in a qualitative 
manner. Finally, an approximate relation between the frequencies observed 
for a certain symmetry type and those calculated by means of force constants 
derived from some other symmetry type is given as an aid in assigning vibrational 
spectra. 


Introduction 


A change in mass, in geometry, or a change in the constants of the potential 
function of a molecule is usually accompanied by a lowering of the molecular 
symmetry. This complicates the direct calculation of the frequency resulting 
from the above changes by making it necessary to solve determinantal 
equations of higher order than for the unchanged molecule. Since calculation 
of the new frequency is so laborious, it is useful to have a simpler method of 
estimating the effect on the frequency of the above mentioned changes. Such 
a method was developed by Lord Rayleigh (4) who devised a perturbation 
theory of vibrating systems to cover such problems. 


If the kinetic and potential energies are expressed in terms of the normal 
co-ordinates @¢, viz: 
2T ad," + Arh.” + asbs” + 
2V = apy + Cobo” + C393? +, 
Lord Rayleigh showed that for changes in a; and ¢; corresponding to da; and dc, 
dd a da; dc 


Mi a, Cy z 
where \y = 42° »? = 42° c 5° for in cm. We may note here that 
dp ,4dx 
ae aD ee, 
D r 


A good many years later in a discussion on the frequency shift due to isotopic 
substitution Teller (7) arrived at the result 
Ox _ LadT; 


Ae ty 
where dT; is the change in kinetic energy associated with substitution of 
atom 7 in the k“ vibration, and 7, is the kinetic energy during the vibration. 


(2) 


1 Manuscript received in original form September 12, 1950, and, as revised, November 6, 1950. 
Contribution from the Division of Chemistry, National Research Council, Ottawa, Canada. 
Issued as N.R.C. No. 2346. 
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rated ‘ ‘ a Bhs Am .. 
By considering the effect of isotopic substitution (when — is small) as a 
m 


perturbation of the kinetic energy Wilson (8) also derived Equation (2). Then 
by treating the problem as a perturbation of the potential energy Wilson (10, 
11) derived an explicit relation between the shift in frequency Av, the co- 
efficients of the perturbing potential, and the coefficients of transformation to 
the normal co-ordinates. In particular, it was shown (9) for Am/m small, 
substitution of one isotope B‘ for B in AB, shifts the totally symmetrical 
vibration by an amount 


=. (3) 


This result had been obtained by Salant and Rosenthal (5, 6) by direct calcu- 
lation of a complex and tedious nature. 


The normal co-ordinates of the system are required to calculate frequency 
shifts by the above mentioned methods of Rayleigh and Wilson. Since it is 
not always convenient to obtain these, the method outlined below is given for 
any set of convenient co-ordinates. The relation (2) of Teller will be derived 
as a special case. 


Change in Frequency Due to Change in Mass (Isotopic Shift) 


Let g: be the displacement of the i" nucleus from the equilibrium position, 
then 


2T = aigr + dogo? + Zaigige +++ (4) 
2V = digi? + boqgo? + Zhigigg +++ . (5) 

The solutions for simple harmonic motion are of the form 
qi= Aix CosV Axt , (6) 


where A;, is the amplitude associated with the i” particle during the k* 
vibration. 


The Lagrange equations of motion are then: 
(b1 — a1 Ax) Are + (bi2 — Gi2 Ax) Azk + + =O (7) 
(bi2 — Giz Ax) Ark + (02 — Ge Ax) Aone + + =O, etc. (8) 
Multiplying (7) by A; and (8) by Asx, etc., and adding it is found that 
nM by Aix? + be Aox? + 2bi2 Aix Aor + + + 
@ Aix” + a2 Ao? + 2a Aix Aoe + + + 
es V(A ix A jx) as, 


= T(Aa Aj) 2 (say) . (10) 


It is possible to expand \ as a function of the amplitudes, b’s and a’s:— 


MA ik + AA ix, a: + Adi, by + Ads) — MA ix, Gi, bi) = Eada ( a) 


ik 


Or Orv, aaa( a) 
Aa; | — Ab; | — ——— J 11 
+ 2Aa (>) +2 (>) +2 2! \OAi? sai 





Ae (9) 
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In all problems which are considered here Aa; and Ad, are small so that higher 
powers may be neglected. In some cases however there are large changes of 
the amplitudes to be considered; as a result the second order term has been 
included in the expansion. Now 











: v 1 ‘ ‘ 
Om _1 ( a Ax ns ) = — X the 7 Lagrange equation 
0A ik T 0A ik 0A ik T 


= 0 
OAK . F : 
and —— = b, — a,x. It is apparent, that for normal co-ordinates, ) is 
ike 

independent of the values of A;,;,* and for generalized co-ordinates is inde- 
pendent of A,, to the first order of small quantities, and for large changes in 
Ax depends on the value of b; — a, Ay. For the k™ vibration the following 
approximations are valid :— 

T~ aR Ake 

vT~ b;. Axe 
so that 6; — ax 4x ~ 0 from (10). Thus, for generalized co-ordinates \ is very 
insensitive to changes in amplitudes. 

Equation (11) now becomes 

A => (2+) Aa; +> 2s) Ab; . (11a) 
a; i 0b; 
For small changes in a; and 6;, substituting for the differentials evaluated from 
(10) into (lla) gives 


dd _ _ Xdrs , Lavy 





r Tk Uk 
oo ee yee (12) 
rs Vi 


This equation reduces to the one (Equation (2)) derived by Teller and Wilson 
when there is no change in the constants of the potential function, i.e., isotopic 
substitution. Equation (12) may be used also to calculate the shift AX/A 
when the energy is expressed in generalized co-ordinates. 


To estimate isotopic shifts it is convenient to adopt the displacement co- 
ordinates 7, r2, 73, etc., and express the kinetic energy in the form: 
2T = nym + nome? + nymsr3?2 +... eee , 
where there are m, M2, 13........ equivalent atoms of mass m, m2, m3...... ‘ 
etc. The potential energy may be expressed then as 
2V = ary? + Core” + 2cwrire + . 
For simple harmonic motion r; = B,, cos (Ax)2t. From Equation (10) 


Ag = 04/Te = Ve/mymM By? + nym2Bo,? + nym3B3.° + , (13) 
where Tk = nym,B,,? + nym2Bo;? + nym3B3) + , (14) 
so d ty, = 1Byx2dm, + n2Bo2dm2 + . (15) 


* Rayleigh’s (op. cit. p. 110) discussion of the stationary property of normal modes suggested 
the above treatment of this problem. 








ARS nm 





ca nti aba ins Mel oan 
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These expressions are valid when all m, m2, m3, etc., atoms have. been substi- 
tuted by the isotopes m', m2‘, m;', etc. If Ni atoms of m, have been substi- 
tuted and N2 atoms of mz and N3 atoms of m3, etc., then 


dt, = NiBij2dm, + N2Boi*dmz + N3B3i2dm;3 + . (16) 
Actual calculation of AX/) falls into the following classes:— 


(i) For 7, a function of m, only, i.e., a stretching vibration, 


(ii) For rt, a function of m, and mz only, i.e., stretching and deformation 
vibrations, 


(iii) When 7; is a function of all the masses, i.e., 1, m2, m3, etc. 


(i) When 7; 1s a function of m, only 


dr = NB, .2dm,/n,By 4.2m, => M ? dm > 
T ny my, 
30 dad _ _ Mdm (17) 
r ny my, 
for as | 
- om 


This is the result obtained by Wilson for the effect of isotopic substitution on 
the totally symmetrical vibration of AB,. 


Some examples and applications of Equation (17) are given in Tables I-VII 
inclusive for CCl, BCl;, NCCI, CsH4, C2H2, HCN, CeHs. The deuterated 
benzenes, ethylenes, and acetylene may be considered as examples since the 
CH group is considered as a whole and thus Am/m is small for substitution of 
a deuterium atom. 


TABLE I 


ISOTOPE EFFECT IN CCl, 


























CCi35 Cci3sCis? | CCis5Cls? | CCI#5C137 | CCl3? 

(Obs)! 461.5 458.4 455.1 

v 1(A1) (Calc)? 461.5 458.2 455.0 451.8 448.4 
(Calc)* 461.5 458.2 455.0 451.8 448.4 
(Obs)! 219.0 216.8 

v 23(E) (Calc)? 219.0 217.4 215.8 214.2 212.6 
(Calc)* 219.0 217.5 216.0 214.5 213.0 
(Obs)! 315.5 312.6 

Vase( Fe) (Calc)3 315.5 313.9 312.3 310.9 308.9 
(Calc)* 315.5 314.4 313.3 312.2 311.1 

| 








1 A. Langseth, Z. Physik, 72: 350. 1931. 

2 Calculated from Equation (17). 

3 Calculated from Equation (18). 

4 Calculated by C. K. Wu and G. B. B. M. Sutherland (J. Chem. Phys. 6:111. 1936.) 
from the equations of Rosenthal, Urey, and Bradley. 
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TABLE II 
ISOTOPE EFFECT IN BCI; 
“| | | 
| B**CR* | Bt °Cl35 B*1Cis*C*" | B*4Cl**Cis* BteCis*Ci** 
“i(A}) (Obs)! | 473.6 473.6 | 469.3 | 464.0 469.3 
(Calc) | 473.6 473.6 | 469.1 | 464.6 469.1 
v(E)| (Obst | 936.7 | 971.0 | | 
(Calc)? | 934.8 | 971.0 | | 
1 B. Bak, J. Jorgensen and Co., Copenhagen (1943). 
2 Calculated from Equation (18). 
3 Calculated from Equation (17). 
TABLE III 
ISOTOPE EFFECT IN NCCI 
| | N¥#CRC] N#CRCI37 | N#CBC]I N#¥CBCIS7 
v». | (Obs)! | 2206 .6 | 2171.0 
|} (Calc)? 2206 .6 2206 .6 2157.0 2157.0 
V2 | (Obs)! | 728.9 7200.8 | 
| (Calc)? 728.9 723.6 706.3 701.3 
1 B. Bak, J. Jorgensen and Co., Copenhagen (1943). 
2 Calculated from Equation (18). 
TABLE IV 
ISOTOPE EFFECT IN ETHYLENE 
| CGH. | GHD | CHD: | GHD, | CD, 
|- —__—__—. — —— 
(Obs)! 1623 1599 j 15724 1545 1515 
Yi | (Calc)? | 1623 1594 1565 1536 1507 
| (Calc) 1623 1596 | 1569 1542 1515 


| 





- 2 mo = 


G. Herzberg, ‘‘Infrared and Raman Spectra,”” D. Van Nostrand Co., N.Y. (1945). 
Calculated from Equation (18) taking CH: as a group. 
Linear interpolation between C2H, and C2D,. 

Average of cis, trans, and asymmetric frequencies. 


TABLE V 


ISOTOPE EFFECT IN ACETYLENE 








| 


| C.H: 








(Obs)! 
(Calc)? 
| (Calc)? 


| 1974 





C:HD C:D2 
1851 1762 
1868 1762 
1933 1882 








eonwrn 


G. Herzberg, ‘‘Infrared and Raman Spectra,’ D. Van Nostrand Co., N.Y. (1945). 
Linear interpolation between C2H2 and C2Dz. 
Calculated from Equation (18) taking CH as a group. 
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TABLE VI 


ISOTOPE EFFECT IN HYDROGEN CYANIDE 











| HCYN" | HC8N" | DCBN" DCYN" 
(Obs)! 2095 =| «2063 | 1926 

V. | 
"| Cale)? 2095 | 2048 2003 2048 











1 B. Bak, J. Jorgensen and Co., Copenhagen (1943). 
2 Calculated from Equation (18) and taking CH as a group. 


TABLE VII 


ISOTOPE EFFECT IN BENZENE 






































C.He C.H;D C.H.D2 C.H;D;3 CsH2D, C.sHD; C.Ds € 
om p|sym vic unsym| 0 m p 

(Obs)! 992.5 | 980.9 | 975.0 962.6 957.9 | 950.9 | 945.2 
11(Aig) (Calc)? | 992.5 | 984.6 | 976.7 968.8 960.9 | 953.0 | 945.2 
(Calc) | 992.5} 985.6 | 978.7 971.8 964.9 | 958.0 | 951.1 

(Obs)! | 1010 1008.6 | 1000 994 980 977 962 

Vi2(Bin) | (Calc)? | 1010 1002 994 986 978 971 962 

(Calc)* | 1010 1003 996 989 982 975 968 

(Obs)! 606.4 | 599 597 593 587 585 579 

vs(E,) | (Cale)? | 606.4 | 602 597 593 588 584 579 
| (Calc)4 | 606.4 | 602.4 | 598.2 594.1 589.9 | 585.7 | 581.5 





1 A. Langseth and R. C. Lord, Kgl. Danske Videnskab. Selskab., XVI; (6) 1. 1938. 
2 Linear interpolation between CysH, and CsDs. 

3 Calculated from Equation (17) taking CH as a group. 

4 Calculated from Equation (18) taking CH as a group. 


(ii) When 7; is a function of m, and m2, one atom of m, and one of m2 are © 
involved in a stretching vibration since these are the terminal atoms of the 
bond. For a diatomic molecule the amplitude is inversely proportional to 
the mass and for a bond in a polyatomic molecule this relation will be approxi- 
mately true, thus 


dr ae M . dm, . M2 
T nN, my, m, + Me 

and so dd ~— M . dm . a E (18) 
r ny my, m, + M2 


For deformation vibrations the masses involved are 
2m, cos? ¢ and my, (apex atom) 


so again the actual masses involved are m,; and mz since for most molecules 
2 cos? @ ~ 1. 
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For bending vibrations then 
d i 9 
Bi eee. a. (19) 
r my mM mM, + Me 


where mz is the apex atom and isotopic substitution is made for m, only. 


Examples of this may be seen in the tables (I-VII). 

(iii) For +, a function of m,, m2, ms, etc., large changes in m; will have a small 
effect on A; when the vibration considered is the one not associated with mass 
m;. From (13) we see that 

const. 





At™~ : 
const. + const. X m3 


Thus 4; ~ a — bm;, (20) 
where a, b are constants, gives the effect on A, of substitution for m3 which is 
not a particle primarily associated with the k“ vibration. 


The figure shows the dependence of the (frequency)? of the nonstretching 




















P . 4 " . 
hydrogen vibrations on the mass = in the compounds XCH2Y. The 
data are taken from Kohlrausch’s book (3). It has been assumed that the 

2.5 
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Fic. 1. Effect of substitution for X and Y in the series X:CH2-Y on the nonstretching 
hydrogen vibrations. 
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potential function for the molecules XCH:2Y has remained unaltered for X, Y 
substitution. From the figure it is seen that the assumptions have been 
justified for the Cl, Br, and I substitution but that for fluorine containing 
compounds a straight line of different slope is obtained. 


Another case of isotopic substitution to be considered is that of deuterium 
for hydrogen, 


H vE D vE 
then Ay = “e and A; = + : 
Tk TR 


but since A is quite insensitive to amplitude changes, 
ee ~ meh - (BY . oN 
a  S é Bix Mpmy : 


Mp 
Cw ao 





(21) 


my 
When T is a function of m only 
gs Mp 
SP ae 


Equation (21) is illustrated in Table VIII where the shift in frequency for 
substitution of one D remains practically constant when additional D’s are 


added. 


TABLE VIII 


LARGE ISOTOPE EFFECT, v;(B,)* 











Molecule Observed Average Deviation from mean 

C.He 3048 

C.H;D 2270 2270 6 
0 2281 

CsH4D2 m 2284 2281 - «5 
p 2279 
sym 2272 

C;3;H;D3; vic 2288 2281 5 
asym 2284 
0 2290 

C.sH2D, m 2272 2276 0 
‘p 2265 

CeHD; | 2268 2268 8 

CDs 2266 2266 10 

2276 4 














* Data and assignment from A. Langseth and R. C. Lord, Kgl. Danske Videnskab. Selskab, 
XVI: (6) 1. (1988). 
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Change in Frequency Due to Configurational Changes 


To compare the frequencies of corresponding modes of vibration for the 
totally symmetrical vibrations of cis and trans dichloroethylene, Equation (10) 
may be of considerable use. The potential function for cis may be written 
2V = 2KcwARcu + 2KecARcci + KecARcé + WuccAPucé + 2Weice 
A@cicc and the kinetic energy is 

2T = 2my(iu” + Ju”) + Wmer(Hci + Her) + Ae(Hc* + Sc”) 
The kinetic energy for the trans molecule has the same expression as that for 
the cis molecule. For the totally symmetrical vibrations the potential energy 
for the trans molecule has the same form as that for the cis but includes a term 
(absent in cis) which is the contribution to the potential energy due to the 
relative motion of the carbon atoms*. If one may assume that interactions 
between nonbonded atoms are small enough to neglect, and that the amplitudes 
of the atoms and force constants in cis and trans molecules are approximately 
equal then 
Nrrans _ Veis + Contribution due to relative CC motion 








Neis Veis 
> 1 (22) 
since the contribution due to the CC motion is +'’*. The frequency of the 
totally symmetrical mode of vibration in the trams molecule is thus greater 
than the frequency in the corresponding totally symmetrical mode for the cis 
molecule. See Table [X. 


TABLE IX 


CHANGE IN FREQUENCY DUE TO CONFIGURATION* 

















| | 
cis trans C2 | Con 
C:H:2Cl. | C2H2Cle C2H.Cl. C.2H,Cl. 
3077 3071 2844 2965 
1587 1576 1445 1437 
A; 1179 Ag 1270 1207 1298 Ay 
711 844 943 990 
173 349 654 748 
265 303 
3072 3080 
B, 1294 B, 1200 
818 817 
571 265 














* Data from H. J. Bernstein and D. A. Ramsay (J. Chem. Phys. 17: 556. 1949.) for she 
dichloroethylenes and H. J. Bernstein (unpublished results) for the dichloroethanes. 


For the antisymmetrical planar vibrations of cis and trans dichloroethylene, 
there is relative CC motion for cis but none for trans, 
trans ps Vtrans 
Neis = Virans + contribution due to relative CC motion 
<1 (23) 


* See F. Lechner, Sitzber. Akad. Wiss. Wien (IIa), 141: 633. 1932, for example. 


thus 











ee es 
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See Table IX for examples. In the case of deformation vibrations where the 
contribution due to relative CC motion is felt most one would expect these 
relations to be valid. For other vibrations however the assumptions may 
break down sufficiently to upset this order (see Table IX for examples). 


A Relation Between the Frequencies Observed for a Particular 
Symmetry Type and the Frequencies Calculated for this 
Type from Force Constants Derived from Some Other Type 


Generally the only symmetry type to which observed frequencies (obtained 
from Raman and infrared spectra) can be assigned unambiguously is the 
totally symmetrical type. The depolarization data from the Raman spectrum 
make this assignment possible. For a molecule of rather high symmetry then, 
there remain usually several more symmetry types to which assignments cannot 
be made from considerations of the spectra only. Usually it is necessary to 
invoke calculations of frequencies based on the force constants calculated 
from the observed totally symmetrical type frequencies. Since there are the 
same number of frequencies as principal force constants in any symmetry type, 
it is necessary to use the simple potential function (i.e., no cross-terms). 
These force constants will in general be first approximations to the true force 
constants which would be obtained if the complete quadratic potential function 
(using all interaction constants) were used. So that using these constants to 
calculate the frequencies belonging to a different symmetry type will at best 
be only a first approximation. Although the numerical agreement between 
the observed and calculated frequencies under these conditions cannot be 
expected to be good, yet there is a very definite relation between the observed 
and calculated frequencies which is of great help in making final assignments. 


From the general relation (10, 11) 
Il = GF, (24) 
where TE = De de De. es sks Ni, 
G = the G matrix which is a function of masses and geometry only, . 
F = the F matrix which isa function of force constants only, we see that 


dil dG, dF 
. 6° FF 
_dF (25) 
F 


-for no mass changes or changes in geometry. 


In the case outlined above where only the principal force constants Ky, 
Ti es K; are used 


F= K,K-2K;3 oeeeee Kes 
Thus UF a1 OF aK, + SF aK, + + \ 
F F \oeK, 0Ke 
dK, , dK» 
=— +++ 





Ky Ke 
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Hence ail = a , drs 4 = dK, 4 dK: 

II Ai Ae Ky, Ke 
If the force constants K, and Kz are calculated from the frequencies observed 
for the totally symmetrical type which has two modes in it, say, then these 
constants used in the secular equation of a different type possessing two modes 
will yield calculated frequencies considerably different from those observed. 
The observed frequencies would correspond to the true force constants which 
may be represented by K; + AK; and Kz + AK». Equation (26) relates the 
per cent difference between the observed and calculated frequencies in terms 
of the per cent difference between the true force constants and those used. 
Since the true force constants are obtained when the interaction constants are 
introduced in the potential function it is evident that upon introduction of a 
small interaction constant the value of dII will be zero. This is readily seen 
by evaluating Ai, 2, from Equations (7), (8) for a symmetry type in which 
there are two modes only, viz.: 


ai (did, + diaz, — 2by2d12) X 4 bib, — bis* oi: (27) 


a\d2 — a" Q\d2 — a2" 





+ . (26) 





Thus introduction of an interaction constant by, say, will appear to the second 
power in the product of the frequencies. To the first order of small quantities 
dil = 0 and Equation (26) becomes 


a Fn dK, 
Ai de Kyi Ke 

The calculated and observed frequencies for this new type will be related in a 

manner such that the per cent increase in one frequency is equal to the per cent 

decrease in the other frequency or vice versa. Also the true force constants 

are related to those used by the same sort of relation 


dK, =. dKe 








Bk 
The problem of three vibrations of a symmetry type for which 


0 = x + ddo dx3 = dK, dKz dK; 





de A3 Ky Ko K; 
reduces to the two aed case if one of the calculated frequencies agrees 
very closely with the value observed (i.e., dA = 0). 








This criterion has been applied (1) with success in the case of the vibrational 
assignment of the Ss molecule which has the symmetry of the D,g point-group. 
From force constants derived from the totally symmetrical class, Bhagavantam 
and Venkatarayudu (2) calculated the frequencies of the type E vibrations to 
be 452, 243, and 122 cm=! Lines at 434, 241, 185, and 152 cm7'! are observed 
in the Raman spectrum. Since these authors employed the criterion of closest 
arithmetical agreement between observed and calculated frequencies their 
assignment for the type £ vibrations was 434, 241, and 151 cm7! However, 
since there is no ambiguity about the line calculated at 452 cm7! and if as an 
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‘approximation it is considered to be very close to the observed value at 434 cm7" 


the problem reduces to estimating the remaining two frequencies by the method 
outlined above. This criterion indicates that two calculated frequencies 
must be chosen such that 


addy = dd» 





Ai Ae 

Of the calculated frequencies 243, 122 and the observed frequencies 241, 185, 
and 152, only 243 and 185, and 122 and 152 satisfy approximately this relation. 
Thus the assignment indicated by this criterion is 434, 185, 152 with 241 
belonging to some other symmetry type. By using this assignment it was 
possible (1) to interpret all of the 35 Raman and infrared bands of Ss, whereas 
with the previous assignment of Bhagavantam and Venkatarayudu about 
eight bands remain unassigned. 
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